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ABSTRACT 

We present improved black hole masses for 35 active galactic nuclei (AGNs) based 
on a complete and consistent reanalysis of broad emission-line reverberation-mapping 
{Sj , data. From objects with multiple line measurements, we find that the highest precision 

measure of the virial product ctAV 2 /G, where r is the emission-line lag relative to 
t^J- \ continuum variations and AV is the emission-line width, is obtained by using the cross- 

correlation function centroid (as opposed to the cross-correlation function peak) for 
i-^ ■ the time delay and the line dispersion (as opposed to full width half maximum) for 

the line width and by measuring the line width in the variable part of the spectrum. 



o 



Accurate line- width measurement depends critically on avoiding contaminating features, 



in particular the narrow components of the emission lines. We find that the precision (or 
random component of the error) of reverberation-based black hole mass measurements 
is typically around 30%, comparable to the precision attained in measurement of black 
hole masses in quiescent galaxies by gas or stellar dynamical methods. Based on results 
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presented in a companion paper by Onken et al., we provide a zero-point calibration 
for the reverberation-based black hole mass scale by using the relationship between 
black hole mass and host-galaxy bulge velocity dispersion. The scatter around this 
relationship implies that the typical systematic uncertainties in reverberation-based 
black hole masses are smaller than a factor of three. We present a preliminary version 
of a mass-luminosity relationship that is much better defined than any previous attempt. 
Scatter about the mass-luminosity relationship for these AGNs appears to be real and 
could be correlated with either Eddington ratio or object inclination. 

Subject headings: galaxies: active - galaxies: nuclei — galaxies: Seyfert — quasars: 
emission lines 

1. INTRODUCTION 

The evidence that active galactic nuclei (AGNs) are powered by gravitational accretion onto 
supermassive black holes is now quite convincing. Certainly there has not yet been a definitive 
detection of the relativistic effects that would be required for unambiguous identification of a 
singularity, although studies of the Fe Ko emission line in the X-ray spectra of AGNs currently 
affords some promise (e.g., Reynolds & Nowak 2002). Nevertheless it seems to be true that the 
centers of both active and quiescent galaxies host supermassive (greater than ~ 10 6 M ) objects 
that must be so compact that other alternatives are very unlikely. 

Black hole masses are measured in a number of ways. In quiescent galaxies, dynamical mod- 
elling of either stellar kinematics (e.g., van der Marel 1994; van der Marel et al. 1998; Verolme et 
al. 2002; Gebhardt et al. 2003) or gas motions (e.g., Harms et al. 1994; Ford et al. 1994; Macchetto 
et al. 1997) is used to determine central masses. In the case of NGC 4258, a weakly active galaxy, 
proper motions and radial velocities of H2O megamaser spots are used to deduce a high precision 
central mass (Miyoshi et al. 1995; Herrnstein et al. 1999). In Type 1 active galaxies (i.e., those with 
prominent broad emission lines in their UV/optical spectra), reverberation mapping (Blandford &: 
McKee 1982; Peterson 1993) of the broad-line region (BLR) can be used to determine the central 
masses. Reverberation mapping is the only method that does not depend on high angular resolu- 
tion, so it is of special interest as it is thus extendable in principle to both very high and very low 
luminosities and to objects at great distances. Moreover, reverberation studies reveal the existence 
of simple scaling relationships that can be used to anchor secondary methods of mass measurement, 
thus making it possible to provide estimates of the masses of large samples of quasars, including 
even very distant quasars, based on relatively simple spectral measurements (e.g., Vestergaard 2002, 
2004; McLure & Jarvis 2002). 

While reverberation methods in principle can be used to determine the full geometry and 
kinematics of the BLR (e.g., Home et al. 2004), applications to date have been comparatively 
simple. Time delays between continuum and emission-line variations are used to deduce the size of 
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the BLR, or more accurately, the size of the line-emitting region for the particular emission line in 
question. By combining the measured time delay r with the emission-line width AV, a virial mass 
can be obtained, 

% , /ctAV 2 

M= J G , (1) 

where c is the speed of light and G is the gravitational constant. The factor / is of order unity and 
depends on the structure, kinematics, and orientation of the BLR. We will sometimes refer below 
to the "unsealed" virial mass M p = ctAV 2 /G, as the virial product, so that the virial mass is the 
virial product times the scaling factor /. 

It became clear even in the first well-sampled reverberation program on NGC 5548 (Clavel et al. 
1991; Peterson et al. 1991; Dietrich et al. 1993; Maoz et al. 1993) that different emission lines have 
different time-delayed responses, or lags. Lags are shorter for lines that are characteristic of more 
highly ionized gases, i.e., the BLR has a stratified ionization structure. It was already known (e.g., 
Osterbrock & Shuder 1982) that higher ionization lines (e.g., HeiiA4686) are broader than lower 
ionization lines (e.g., H/3 A4861), and it was natural to look for a virial relationship between lag and 
line width, r oc AV~ 2 , which would constitute evidence that gravity dominates the motions of the 
BLR gas and that the black hole mass can therefore be inferred. Early attempts to do this were not 
promising, although Krolik et al. (1991) did note the trend of decreasing time lag with increasing 
line width for the UV lines in NGC 5548. However, upon revisiting the issue, Peterson & Wandel 

(1999) found that there is indeed a virial relationship between lag and line width in the case of NGC 
5548 if the line width is measured in the variable part of the emission line and one avoids (a) lines 
that are strongly blended with other features and (b) lines with lags that are uncertain because of 
potential aliasing effects in the time-series analysis. Similar virial-like relationships between lag and 
line width were subsequently found in other objects (Peterson & Wandel 2000; Onken & Peterson 
2002; Kollatschny 2003). Despite earlier claims that emission-line reverberation yielded masses that 
were too low by a factor of several (Ho 1999), it was subsequently shown (Gebhardt et al. 2000b; 
Ferrarese et al. 2001) that the relationship between AGN reverberation-based black hole masses 
Mbh and their host-galaxy bulge velocity dispersions appears to be consistent with the black 
hole mass/bulge velocity dispersion relationship (hereafter the Mbh~c* relationship; Ferrarese k, 
Merritt 2000; Gebhardt et al. 2000a) that is seen in normal galaxies. Moreover, the relationship 
between black hole mass and host-galaxy bulge luminosity also seems to be the same for both 
quiescent and active galaxies (Wandel 2002; McLure & Dunlop 2002). 

Unfortunately, there is a significant systematic uncertainty (Krolik 2001) in AGN reverberation 
masses embodied in the scaling factor / in eq. (1), which remains unknown. For lack of a better 
estimate, published studies have usually used a nominal value of / = 0.75 for AV taken to be the 
full width half maximum (FWHM) of the emission line, as described in section 6 below. 

Wandel, Peterson, & Malkan (1999; Paper I in this series and hereafter referred to as WPM) 
published a compilation of black hole masses in 17 Seyfert 1 galaxies and two quasars. Kaspi et al. 

(2000) published the results of a large reverberation-mapping campaign that led to mass measure- 
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ments for the central objects in 17 Palomar- Green (Schmidt & Green 1983) quasars and combined 
their results with WPM to obtain relationships between the BLR size and AGN optical luminos- 
ity (the "radius-luminosity relationship") and between the central mass, which we will henceforth 
assume to be a black hole, and the optical luminosity (the "mass-luminosity relationship"). These 
are both of obvious importance: 

1 . The radius- luminosity relationship can be used to deduce the masses of black holes in distant 
quasars by combining the inferred BLR radius with the widths of the emission lines. 

2. The mass-luminosity relationship hence obtained relates directly to current accretion rates 
and radiative efficiencies. The mass itself provides a strong constraint on the black hole 
growth history. 

All of the important relationships mentioned here — the time-lag/line-width virial relation- 
ship, the AGN Mbh-c* relationship, the radius-luminosity relationship, and the mass-luminosity 
relationship — show considerable scatter. Moreover, the reverberation database is very inhomoge- 
neous, and the data have not always been analyzed in a uniform way; this is particularly true in 
the case of many of the earliest results. We suspected that more than a decade of experience in de- 
veloping reverberation mapping techniques and error analysis merited reanalysis of the earlier data 
and, that in at least some cases, improved calibration for spectra would lead to improved results. 
This suspicion was borne out in the case of NGC 3783 (Onken & Peterson 2002); our reanalysis of 
the combined UV and optical data led to a determination of the virial mass of the central object 
that was an order of magnitude more precise than that quoted by WPM, which was based only on 
the original optical spectra and optical continuum and H/3 light curves. Fundamentally, a complete 
reanalysis of the body of reverberation database is warranted by relatively recent (a) improvements 
in cross-correlation error analysis, (b) recognition of the importance of measuring line widths in 
the variable part of the spectrum, and (c) recognition that emission-line time lags can vary over 
time scales longer than the reverberation time scale due to changes in the mean luminosity of the 
object (Peterson et al. 2002). 

We thus decided to undertake a massive reanalysis of all of the reverberation mapping data 
readily available to us, for the express purpose of improving AGN black hole mass determinations. 
We distinguish in the usual way (e.g., Bevington 1969) between the accuracy of these masses (i.e., 
how close they are to the true values), which depends on how well we can account for systematics, 
and their precision (i.e., how exactly we measure the virial products), which depends primarily on 
"random errors" associated with measuring line widths and time lags. Thus, this investigation has 
two parts, with different goals: 

1. In order to improve the precision of AGN black hole mass measurements, we are reanalyzing 
all of the readily available reverberation data to determine the best measures of time lag and 
line width for these studies. We do this by assuming that the virial product for all emission 
lines is the same at all times for a particular AGN. We then explore ways of characterizing 
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lags and line widths that yield relationships most consistent with a virial relationship. This 
is the subject of this contribution. 

2. In order to improve the accuracy of AGN black hole mass measurements, we are obtaining 
high-precision measurements of cr* for reverberation-mapped AGNs. We demonstrate that 
there is general consistency between the Mbh - c* relationships between AGNs and quiescent 
galaxies, and then assume that these two relationships have a common zero point, thus 
determining the scale factor /. This is the subject of a companion paper (Onken et al. 2004), 
which we will draw on for the absolute calibration of the black hole mass scale. 

2. DATA 

We have included in this analysis all objects for which we had ready access to the spectra used 
in the original investigations. This consists of data from most of the reverberation-mapping exper- 
iments undertaken to date, including the large samples from International AGN Watch 1 projects 
(Alloin et al. 1994; Peterson 1999), the Lovers of Active Galaxies (LAG) campaign (Robinson 1994), 
the Ohio State monitoring program (Peterson et al. 1998a), and the Wise Observatory /Steward 
Observatory monitoring program (Kaspi et al. 2000). A list of objects analyzed here is given in 
Table 1. Column (1) gives the common name of the object as used in the relevant papers on the 
reverberation results, which are referenced in column (2). Epoch 2000 coordinates are given in 
columns (3) and (4). The redshift z of each object is in column (5), with nominal V magnitude 
and 5-band extinction, based on 100 /im dust maps from Schlegel, Finkbeiner, & Davis (1998), are 
given in columns (6) and (7), respectively. Column (8) gives the standard name of the object in 
the Veron-Cetty Sz Veron (2001) catalog, and column (9) gives other common names by which the 
object is often known in the AGN literature. All entries in this table are from the NASA/IPAC 
Extragalactic Database. 

For this analysis, the fundamental data that we require are (1) continuum and emission-line 
light curves to determine time lags and (2) spectra from which line widths are to be measured. 
We use the published versions of the light curves, except where noted below. For the purpose of 
measuring the line widths, we use the same spectra from which the continuum and emission-line 
measurements were made; in some cases, notably the ground-based component of the International 
AGN Watch monitoring programs on Fairall 9, NGC 3783, NGC 4051, NGC 4151, NGC 5548, 3C 
390.3, and NGC 7469, we restricted consideration to the single largest homogeneous data subsets, 
i.e., those data which are most similar in terms of resolution and quality, often from a single source. 



x Data obtained as part of International AGN Watch projects are available at http://www.astronomy.ohio- 
state.edu/~agnwatch/. 
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For each set of spectra, we formed a mean spectrum, 



(2) 



where -Fj(A) is the zth spectrum of the N spectra that comprise the database. We also define a 
root-mean-square (rms) spectrum 



We form the mean and rms spectra only from the most homogeneous subsets of the database, taking 
care to ensure that variations observed in the homogeneous subset are consistent with those observed 
in the entire data set. Strict homogeneity of the data, particularly in terms of spectral resolution 
and spectrograph entrance aperture, is necessary to avoid introduction of spurious features in the 
mean and rms spectra. 

Line widths can be measured in either the mean or rms spectrum; the advantage of using 
the rms spectrum is that constant components of the spectrum, or those that vary on timescales 
much longer than the duration of the experiment, vanish, thus largely obviating the problem of 
deblending lines. The corresponding disadvantage, however, is that the rms spectra are generally 
much noisier than the mean spectra as the amplitude of variability is usually fairly small for these 
AGNs. The most compelling reason to use the rms spectra is that then we are measuring the parts 
of the emission lines that are actually varying. We thus have a strong prejudice towards using the 
rms spectra, and will attempt to justify this choice below. 

For several of the galaxies listed in Table 1, there are multiple data sets available, sometimes 
from the same source and sometimes from different sources. We analyze each individual set as 
an independent time series. In the case of some of the brighter Seyfert 1 galaxies which have 
relatively short H/3 response times, multiple- year campaigns were broken down into individual 
subsets covering single observing seasons, thus yielding multiple, independent measurements of the 
line widths, lags, and virial masses. This is desirable not only from a statistical point of view, but 
because it is now clear (Peterson et al. 2002) that both lags and line widths can vary from one 
observing season to the next as the mean continuum luminosity slowly varies. 



In Table 2, we identify each individual data set for which time-series analysis was carried 
out. Column (1) gives the common name of the object, and column (2) gives the reference for the 
original data. Individual emission lines are identified in column (3) and the spectral resolution of 
the data (see below) appears in column (4). Column (5) gives the range in Julian dates spanned 





3. 



LINE WIDTHS 
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by the spectra 2 . Individual emission lines were isolated by interpolating a linear continuum (in 
units of ergs s _1 cm~ 2 A -1 ) between continuum windows on either side of the line (columns 6 and 
7) underneath the line, whose limits are given in column (8); all wavelengths in Table 2 are in the 
observer's reference frame. 



3.1. Measures of Line Width 

Given an emission- line profile P(X) (i.e., flux per unit wavelength above a continuum interpo- 
lated underneath the line), we parameterized the line width in two separate ways: 



FWHM. How this quantity is measured depends on whether the line is single or double peaked. 
In the case of a single-peaked line, we identify the line peak -P(A) max . We then start at the 
short-wavelength limit of the line (column 8 of Table 2) and search for the Ai such that -P(Ai) = 
0.5P(A) 

max- We then repeat the search starting from the line peak and moving to shorter wave- 
lengths to find A2 such that -P(A2) = 0.5i- > (A) max . The mean of these two wavelengths is taken to 
be the wavelength at half-maximum flux on the short-wavelength side of the profile. An identical 
procedure is used to identify the half-maximum point on the long-wavelength side of the line, and 
the difference between these is taken to be the FWHM. For a double-peaked line, we define a 
short-wavelength peak P(A) max , s hort and a long-wavelength peak -P(A) maXj i ong . We then follow pro- 
cedures similar to those above: we define Ai and A2 relative to the short-wavelength peak only, and 
compute their mean. A similar calculation is done on the long-wavelength side, this time relative 
to -P(A) maXi i ong , and the FWHM is taken to be the separation between the calculated means of Ai 
and A2 on either side of the line. This is illustrated in Fig. 1. 



Line Dispersion. The first moment of the line profile is 

A = J XP(X)dX/ J P(X)dX. (4) 

We use the second moment of the profile to define the variance or mean square dispersion 

^ine(A) = (A 2 ) " A o = (/ X 2 P(X)dX/ J P(X)dX^j - A 2 . (5) 

The square root of this equation is the line dispersion aii nc or root-mean square (rms) width of the 
line. 

Both measures have intrinsic strengths and weaknesses: FWHM is trivial to measure, except 
in the case of multiple-peaked lines or noisy data, and can even be accurately estimated graphically. 



2 It should be noted that in many cases, isolated points at the beginning or end of the original time series may be 
excluded from our analysis. 
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Compared to crime ) it is also less sensitive to blending with other lines and to the contribution from 
extended line wings. On the other hand, <7ii ne is well defined for arbitrary line profiles, less sensitive 
to the presence of even fairly strong narrow-line components, and as we shall see throughout this 
analysis, more accurate for low-contrast lines and the relative uncertainties are much lower than 
for FWHM. On the other hand, oii ne is in some cases also problematic; an nc — ► oo for a Lorentzian 
profile, for example (though in practice, the wings of any reasonable line profile become lost in 
the noise). Fromerth & Melia (2000) point out many of the advantages of <7ii ne relative to FWHM 
and that show a virial relationship between lag and line width is also found using a"ii ne rather than 
FWHM to characterize the line width. 

It is worth noting at this point that there is a simple relationship between these two quantities 
for a given line profile. For a Gaussian, FWHM/crn ne = 2(2 In 2) 1 / 2 ~ 2.355, and for a rectangular 
profile (produced by emission-line clouds in circular Keplerian orbits of fixed radius and random 
inclination), FWHM/ai ine = 2(a) 1 / 2 » 3.464. 

3.2. Resolution Correction 

Since some of the emission lines widths are actually rather narrow, we need to correct each line- 
width measurement for the finite resolution of the spectrograph with which the data were obtained. 
We will assume that the observed line widths AA Q b s can be written in terms of the intrinsic line 
widths AAtruc and the spectrograph resolution AA rcs as 

AA 2 bs « AA t 2 mc + AA r 2 cs . (6) 

Application of this equation to obtain AAtme requires knowing the resolution at which the obser- 
vations were made. In order to determine this for the optical data used here, we relied on accurate, 
high-resolution measurements of the width of the [O m] A5007 line in many of the AGNs discussed 
here by Whittle (1992). Whittle's FWHM measurements for AGNs in this study are given in Table 
3, in the rest frame of each galaxy. In order to determine the resolution of the data used in this 
study, we transformed the values in Table 3 back to the observed frame and to wavelengths units 
and assumed this to be AAtme- We then took our measurements of FWHM([0 in] A5007) as AA Q b s 
and solved for AA res , the FWHM resolution of the data. These are the values given in column (4) 
of Table 2. We recover AA true for the broad lines by application of eq. (6). 

The absence of isolated narrow lines in the UV spectra of AGNs precluded using this method 
for UV spectra. Instead, we assume a spectral resolution of 6 A for the International Ultraviolet 
Explorer (WE) SWP camera and a resolution of 1.9 A for Hubble Space Telescope (HST) FOS 
spectra (e.g., Korista et al. 1995). 
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3.3. Narrow-Line Contamination 

As noted above, in principle constant line components should not appear in rms spectra. In 
practice, we find that residual narrow-line features, generally weak but sometimes quite strong, 
appear in our rms spectra. These narrow-line residuals appear when the data are less than ide- 
ally homogeneous, in particular, when the line-spread function is not the same for each spectrum. 
The Wise Observatory spectra are particularly prone to this because of the method used to effect 
a photometric calibration of the spectra; these data are obtained through a long, wide slit that 
also contains a nearby field star that is used for relative photometric calibration. An unfortunate 
side effect of this highly accurate photometric calibration method is that the line-spread function 
is not well-controlled as the target can migrate small amounts within the slit in the dispersion 
direction; while this has only a small effect on the measured broad-line widths, the narrow-line 
profiles are strongly affected and fairly strong narrow-line residuals can result in the rms spec- 
trum. Again as we will see below, accurate narrow-line removal is critical for accurate line-width 
measurement. Therefore, in cases where the rms spectra show significant residual narrow compo- 
nents of H/3 and/or [O in] AA4959, 5007, we have subtracted these components from each spectrum 
before combining them into mean and rms spectra. In such cases, we used [O m] A5007 as a 
template profile and shifted and flux-scaled this profile to obtain a suitable model of the narrow 
lines, which we then subtracted from each spectrum prior to forming the mean and rms spec- 
tra. The [Oiii] A5007/[Oin] A4959 flux ratio is fixed at a value ~ 3 by atomic physics, but the 
H/3 ri arrow/[0 in] A5007 flux ratio is different for each galaxy. For galaxies in which narrow-line re- 
moval was undertaken for even some of the data, the adopted narrow-line fluxes are shown in Table 
4. Most of these values are from Marziani et al. (2003), although in few cases we used our own 
determinations. 

Unfortunately, decomposition of the narrow lines from the spectra is much more difficult in 
the Ha and ITy regions of the optical spectrum, and was therefore not attempted. Cases in which 
this might present a problem are noted below. Finally, we note that none of the UV data from IUE 
or HST show narrow-line residuals, consistent with little or no narrow-line contribution to the UV 
emission lines in quasars (cf. Wills et al. 1993). 

3.4. Line Width Uncertainties 

To determine FWHM and <7ii ne and their associated uncertainties, we employ a bootstrap 
method similar to that which we use for the time-series analysis. A given data set contains N 
spectra. For a single bootstrap realization, we make N random selections from this group, without 
regard to whether or not a particular spectrum has been previously selected. From these N spectra 
we form a mean spectrum and an rms spectrum (eqs. 2 and 3), and fits to the underlying continuum 
and the line measurements are performed on these spectra. Multiple bootstrap realizations allow 
us to build up a distribution of line- width measurements from the random sets of N spectra. From 
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these we can compute a mean value and standard deviation for the width of each emission line, 
and these are the values that we will use in this analysis. 

This method of determining the line widths and associated errors is different than what we have 
done previously. Line widths and uncertainties presented by WPM were determined less rigorously, 
by comparing the measurements of FWHM obtained by using the "highest plausible" and "lowest 
plausible" underlying continua. Figure 2 shows a comparison between FWHM values and associated 
uncertainties from this study and those reported by WPM; note that for this particular comparison 
only, we did not adjust our measurements for spectral resolution in order to effect a more meaningful 
comparison with WPM. In general, the measurements and errors are both in good agreement. The 
uncertainties we find here are on average about 12% lower than those quoted by WPM. 

The results on two quasars in Table 1, PG 1351+640 and PG 1704+608, were deemed to be 
too poor to retain in this analysis; the emission- line variability was simply too weak to produce 
reliably measurable emission lines in the rms spectra. We note in particular that the problem with 
PG 1704+608 has already been discussed in the literature (Boroson 2003); the H/3 line in the rms 
spectrum is not in fact variable broad-line emission, but merely residual narrow-line emission. 

4. TIME-SERIES ANALYSIS 

The methodology we employ for measurement of time lags and their associated errors is the 
interpolation cross-correlation method, essentially as described by White & Peterson (1994) and 
by Peterson et al. (1998b, hereafter P98b), but with some modifications that are described in 
the Appendix. A complete tutorial on our cross-correlation methodology is provided by Peterson 
(2001). The estimates of the uncertainties are based on a model-independent Monte Carlo method 
in which a single realization yields a cross-correlation function (CCF) whose centroid r ccnt , peak 
value r max , and peak location r pea k are measured. As discussed in the Appendix, we compute r ccn t 
using only the points at values r > 0.8r max , where r max is the peak value of the CCF. A large 
number of independent realizations is used to build up a "cross-correlation centroid distribution 
(CCCD)" and a "cross-correlation peak distribution (CCPD)" (cf. Maoz & Netzer 1989). We take 
T C ent and Tp Ca k to be the means of these distributions. The CCCDs and CCPDs are generally non- 
Gaussian, so we define upper and lower uncertainties separately such that 15.87% of the realizations 
yield values larger than the mean plus the upper error and that 15.87% of the realizations yield 
values smaller than the mean minus the lower error (i.e., the errors are ±1<7 errors if the distribution 
is Gaussian). 

We carried out a cross-correlation analysis for each dataset, as summarized in Table 5. The 
object is listed in column (1) and the emission line and relevant Julian Date range are listed in 
columns (2) and (3), respectively (cf. Table 2). Column (4) gives the peak amplitude of the CCF 
r max; an d columns (5) and (6) give the noise-corrected fractional variation F var (cf. Rodriguez- 
Pascual et al. 1997) of the continuum and line, respectively, during the range of dates given in 
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column (3). Columns (7) and (8) give the CCF centroid T ccn t and CCF peak r pea k, respectively, both 
in the observed frame. Uncertainties in these quantities were estimated by employing the model- 
independent FR/RSS method of P98b, with selected modifications suggested by Welsh (1999), as 
described in the Appendix. The uncertainty associated with r max (column 4) is the rms variation 
in this quantity for the multiple Monte Carlo realizations. Note that entries preceded by colons 
(:) are deemed to be unreliable (see the appropriate notes on the individual objects below), i.e., 
these are cases where there may be systematic uncertainties larger than indicated by the quoted 
uncertainties. 



5. TESTS OF VIRIAL RELATIONSHIPS 

Our next goal is to determine empirically which measures of time delay (r cen t or Tp ea k) an d 
line width (FWHM or cii nc ) provide the most robust estimates of the black hole masses. Specifi- 
cally, we consider which combination of these measures gives us the most consistent or minimum 
variance virial product, ctAV 2 /G, where r is the time delay and AV is the line width. This test 
can be performed on four of the objects in our sample for which long-duration multiwavelength 
spectroscopy allow measurements of a number of different variable emission lines. By far the best 
and most extensive data are those on NGC 5548, and we give these great weight in our analysis. 
On the other hand, we do not give much weight to the results for 3C 390.3 on account of relatively 
large uncertainties in both time lags and line widths. 

Table 6 gives our measurements for the time-lag and line-width parameters. Columns (1) and 
(2) identify the object and data set, in the same order as in Tables 2 and 5. Columns (3) and 

(4) contain r ccnt and r pca k, respectively; these are the values in Table 5, now corrected for time 
dilation 3 by dividing by 1 + z. Line-width measurements were transformed to the rest frame of the 
object and converted to line-of-sight velocities; the values for aii nc and FWHM are given in columns 

(5) and (6), respectively. Again, a colon preceeding an entry indicates that we do not regard the 
entry to be reliable. 

At this point, we make the assumption that the most robust measures of the time delay and 
line width are those that most closely yield the virial relationship AV oc r -1 / 2 . The justification 
for this assumption is simply that a virial relationship between time delay and line width has 
already been established for several objects (Peterson & Wandel 1999, 2000; Onken Sz Peterson 
2002; Kollatschny 2003). We proceed by examining the four cases where multiple measurements of 
the virial product are available, giving the most weight to the results on NGC 5548. 



3 WPM did not apply a time-dilation correction for their objectss as the redshifts are all low. Kaspi et al. (2000) 
made this correction for their own higher-redshift objects and the objects of WPM. 
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5.1. Virial Relationships in Individual Objects 

NGC 5548. Figure 3 shows four plots of the virial relationship for all the lines in NGC 5548, as 
measured in the rms spectra. The optical data from the International AGN Watch (Peterson et al. 
2002 and references therein) are divided into subsets based on single observing seasons separated 
by the several-week gap when NGC 5548 is too close to the Sun to observe. Some experimentation 
revealed that division into shorter subsets led to much larger errors in lag measurements. The four 
panels in Fig. 3 show the possible permutations of the virial product using r ccn t and T pea k for the 
time-delay parameter and uii ne and FWHM as the line-width parameter. Fits to these data and 
the other virial relationships described below are summarized in Table 7, in each case for a best-fit 
slope and for a force-fit to a slope of b = —1/2. Fits were obtained using the orthogonal regression 
program GaussFit 4 (Jefferys, Fitzpatrick, & McArthur 1988), which accounts for errors in both 
parameters. These data show clearly that (a) the virial relationship AV oc r -1 / 2 is robust, i.e., 
precisely how the time delay and line width are measured is not critical, and (b) the least scatter 
in this relationship is obtained by using r ccnt and cri inc to parameterize the relationship. This is 
confirmed by computing the virial product ctAV 2 /G for each measure; the combination of T ccn t 
and uiine has a mean precision (standard deviation divided by the mean) of about 0.032, which is 
lower than for the other pairs of measurements. 

It is important to point out that the fractional errors in FWHM are rather larger than those 
in C7ii ne . A simple consequence of this is that the x 2 statistic can be misleading, as it is larger for 
o"ii nc than for FWHM. In any case, we hasten to point out that the scatter in these relationships is 
sufficiently large that it is clear that simple virial motion is an incomplete description of gas motions 
in the BLR. Figure 4 shows the virial relationship obtained by using the line widths measured from 
the mean spectra, uncorrected for narrow-line contamination. The deleterious effect of the narrow- 
line contribution on the line width measurements is most strongly apparent for H/3, as expected. 
In Fig. 5, we show that correcting the spectra for narrow-line H/3 improves the result, but only 
somewhat. It seems clear that the rms spectrum should be used for these measurements. 

NGC 3783 and NGC 7469. Figures 6 and 7 show the virial relationship for NGC 3783 and 
NGC 7469 respectively, two of the other well-studied AGNs for which multiple measurements of the 
emission-line lags are available, though in both cases there are far fewer data than for NGC 5548. 
The measurements for NGC 3783 are in excellent agreement with viral relationship. The results 
for NGC 7469 are in poorer agreement, but the best-fit slope is within 2a of the virial prediction. 
Again our conclusions do not hinge critically on which time-lag and line-width measures we use. 

3C 390.3. Figure 8 shows the case of 3C 390.3, which seems to afford some difficulties. The left- 
hand column shows a plot of line-width measures, <7ii ne in the top panel and FWHM in the bottom 



4 GaussFit is publicly available at ftp://clyde.as.utexas.edu/pub/gaussfit/. 
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panel, versus r cen t- The apparent lack of consistency with a virial relationship arises because: 

1. All the lines in this object are very broad and each of the measured lines are contaminated 
by blending: Lya is blended with NvA1240, CivA1549 is blended with HeilA1640, and 
He II A4686 and H/3 are blended. 

2. The emission-line lags have relatively large uncertainties and span a comparatively limited 
range (the largest and smallest lags differ by a factor of only ~ 1.5, compared to a factor of 
7-14 for the other three objects discussed above). 

In an attempt to circumvent the problem of line blending, we will assume that each line is intrin- 
sically symmetric about its nominal wavelength. Its true width therefore can be better estimated 
by measuring only the unblended half of the line and reflecting it about the line center. The virial 
relationship using these line widths is replotted in the right-hand column of Fig. 8. This gives 
somewhat improved consistency with the virial relationship, especially for FWHM. However, the 
relatively low range in lags and line widths and large errors in lag measurements hardly make this 
object a convincing case for a virial relationship. Given the large systematic uncertainties in the 
line widths on account of line blending and the relatively poor precision of time-delay measurements 
we do not believe that these results are inconsistent with a virial relationship. 

5.2. General Results 

We conclude from this analysis that the most consistent virial product is obtained by using 
T C cnt and criine as the time-lag and line-width measures. The discussion in the Appendix further 
assures us that T ccn t is a good choice for the lag measurement. The virial product computed from 
T C cnt and uiine is thus given for each data set in column (7) of Table 6. 

We also conclude that for the purpose of determining black hole masses, the rms spectrum 
provides the most reliable line width measurement. However, it is also clear (e.g., Fig. 5) that 
the mean spectrum (or perhaps even a single spectrum) , with its much higher signal-to-noise ratio 
than the rms spectrum, can be used with little penalty in accuracy, as long as one can adequately 
account for contamination by other features, notably the narrow- line contribution and blending 
with adjacent features. However, the strength of narrow-line contributions to broad-line spectra 
are often known rather poorly if at all and blending by various features, notably Fe n contamination 
of H/3, is problematic. Further discussion is beyond the scope of the current paper, but will be 
pursued elsewhere. 

For completeness, we show in Figs. 9 and 10 the distribution of the ratios FWHM/crii ne and 
Tp ea k/T ccn t, respectively, for all lines used in this analysis (i.e., highly uncertain values excluded). 
The mean and standard deviations of these distributions are FWHM/(Ji inc = 2.03 ± 0.59 and 
r peak/fccnt = 0.95 ± 0.20. The low mean value of FWHM/cri ine relative to that for a Gaussian 
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(FWHM/(Tii n e = 2.355) means that, on average, these lines have weaker cores and stronger wings 
relative to Gaussians. 



6. BLACK HOLE MASSES 

For many of the objects in this study, we have multiple measurements of the virial products. 
In Table 8, we list for each object in this study the weighted mean virial product < CT ccn t<jf ine / G >, 
based on the entries in column (7) of Table 6, excluding the more uncertain values (those preceeded 
by a colon). 

The masses of the central objects are given by 

M BH = /CTcc ° taL % (7) 

where as noted earlier / depends on the structure, kinematics, and aspect of the BLR. The scaling 
factor / can be determined in a number of ways, the easiest being to assume that AGNs and 
quiescent galaxies follow the same Mbh _ <7* relationship; one can then use quiescent galaxy results 
to normalize the AGN Mbh-c* relationship and provide an absolute mass scale for AGN black 
holes. We carry out this exercise in a companion paper (Onken et al. 2004), in which we find 
(/) = 5.5. Our final black hole masses, based on eq. (7) with an adopted mean value (/) = 5.5, are 
given in column (3) of Table 8. 



6.1. Uncertainties in Black Hole Masses 

As noted earlier, the first goal of this project has been to improve the precision of the virial 
product measurement. We find that the typical precision (i.e., fractional error) of the virial product 
measurement is about 33% for the 35 AGNs for which we are able to estimate black hole masses, or 
~ 26%, excluding NGC 4593 and IC 4329A, for which the reverberation results are notably poor. 
The second goal is to improve the statistical accuracy of the reverberation-based black hole mass 
scale using the normalization of the AGN Mbh _ <7* relationship reported in a companion paper 
(Onken et al. 2004). The scatter around the AGN Mbh~c* relationship is found to be about a 
factor 2.6 - 2.9, depending somewhat on the slope of the quiescent galaxy Mbh~<7* relationship 
(Onken et al. 2004). It is important to keep in mind that this level of accuracy is statistical in 
nature and individual black hole masses may be less accurate. 

It must be kept in mind that there are various systematic difficulties with reverberation results 
that can in principle lead to significant errors in individual black hole masses (e.g., Krolik 2001). 
Reverberation, of course, is not the only method of measuring black hole masses that can fail 
catastrophically under certain conditions: both stellar dynamical and gas dynamical methods can 
also lead to ambiguous or even misleading results (e.g., Verdoes Kleijn et al. 2000; Cappelari et al. 
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2003; Valluri, Merritt, & Emsellem 2004). While we find the limited scatter in the AGN M BH -a* 
relationship reassuring, additional tests remain highly desirable. 

6.2. On Normalization of the AGN AIbh — f* Relationship 

As noted above, we normalize the AGN Mbh~c* relationship to the Mbh _ c* relationship for 
quiescent galaxies by setting / = 5.5 in eq. (7). This represents the first empirical determination 
of the zero point in the AGN black hole mass scale. 

In previous work, the scale factor appeared in a different and more model-dependent way, only 
as an adjustment e to the line-width parameter, i.e., 

a = eFWHM. (8) 

and the virial mass is M = c{t)o 2 /G. For example, Netzer (1990) assumed e = \/3/2, which arises 
from assuming that cr line = FWHM/2 and that the velocity dispersion is isotropic, i.e., a = valine) 
recalling that <7ii ne is the line-of-sight velocity dispersion. This further assumes that the mean time 
delay (r) is independent of aspect or inclination, which is true when the line emission is isotropic 
and unabsorbed and the geometry has spherical or polar symmetry. If we make the equivalent 
assumptions in eq. (7), i.e., that the velocity dispersion is isotropic so a = y/Sa\i ne and that 
(June ~ FWHM/2 (which is on average quite a good assumption, given the results of the previous 
section), then / = 3. Our empirical calibration of the AGN mass scale through the AGN Mbh~ct* 
relationship is thus about a factor of 1.8 (= 5.5/3) higher than the mass scale used in previous 
papers (e.g., WPM and Kaspi et al. 2000). We emphasize again that the previous normalization 
was made only in the absence of observational information or better-justified assumptions; the value 
we give here is the first determination based on observational parameters, in this case as embodied 
in the Mbh _ c* relationship. 

7. COMMENTS ON SELECTED INDIVIDUAL OBJECTS 

PG 0026+129. Because of strong narrow-line residuals in the rms spectrum, we found it neces- 
sary to remove the narrow component of H/? and the [O III] AA4959, 5007 lines from the individual 
spectra. However, the narrow lines in the vicinity of Ha are sufficiently weak that they present no 
problems. 

PG 0052+251. The narrow-line components were removed from the H/3 region, but the Ha 
region suffers from significant narrow-line residuals in the rms spectrum. The FWHM of Ha is 
untrustworthy and probably underestimated for this reason. 
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Fairall 9. In this case, we used the time-binned UV continuum light curve at 1390 A (Rodrfguez- 
Pascual et al. 1997) as the driving light curve for all lines, including H/3. Unfortunately, the 
C iv A 1549 light curve is very noisy, and the cross-correlation result is not trustworthy. It is therefore 
excluded from the mass determination. 

Mrk 590. The H/3 profile in the rms spectrum in the first data set (JD2448090-JD2448323) seems 
anomalously narrow, although the rms spectrum for this data set is significantly noisier than for 
the other rms spectra (note the lower value of F var in Table 5). We have therefore excluded the first 
data set from the mass analysis, since the line width is questionable although the lag measurement 
seems trustworthy. 

Mrk 79. The lag measurements from the fourth data set (JD2449996-JD2450220) are not trust- 
worthy on account of some significant aliasing effects, seen clearly in the double-peaked CCCD 
shown in Fig. 11. We have therefore excluded this data set. 

PG 0804+761. This is one of a handful of objects that are similar to the well-known Seyfert 
galaxy I Zw 1. These "I Zw 1-like" AGNs have relatively narrow lines and very strong optical Fen 
emission; they are among the more extreme members of the narrow-line Seyfert 1 (NLS1) subclass. 
A particular problem these objectss present is that the [O in] AA4959, 5007 lines are heavily con- 
taminated by Fen emission, and in fact, most of the emission that makes up these features can be 
attributed to FeiiA4924 and FeiiA5018 (e.g., Peterson, Meyers, & Capriotti 1984). We are thus 
unable to use [O in] A5007 as a template for removal of the narrow-line contaminants. However, in 
these cases, narrow H/3 is usually too weak to strongly affect the line- width measurements, so we do 
not attempt narrow-line removal for these objects. In this particular quasar, there is a clear residual 
narrow-line Ha component in the rms spectrum, but the line-width measurement is probably not 
strongly affected. On the other hand, a strong [O in] A4363 residual makes H7 highly asymmetric 
in the rms spectrum, and thus the line-width measurements for H7 cannot be trusted. 

PG 0844+349. This is another I Zw 1-like object (see PG 0804+761 above). The time-lag 
measurements are very inconsistent from line-to- line, probably because of inadequate time sampling. 
There are narrow-line residuals in H7 and Ha. The H/3 cross-correlation function is clearly strongly 
affected by correlated errors and is thus rejected. The results on this object are of rather low quality. 

Mrk 110. The H/3 data used here are from Peterson et al. (1998a). HeiiA4686 appears as a 
prominent, broad feature in the rms spectra of this object, as shown in Fig. 12. We therefore con- 
structed a light curve for He 11 from the original data. Unfortunately, the He 11 lags are so short that 
the measurements of r cen t and r pea k cannot be trusted, as they are significantly shorter than the time 
interval between observations. We therefore do not use the He 11 lines in the mass determination, 
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although we point out that they are generally consistent with the H/3 results. Kollatschny (2003) 
has also studied the variability of this AGN, and finds that the results for Ha, H/3, He I A5876, and 
He II A4686 are consistent with a single virial mass. As discussed in more detail by Onken et al. 
(2004), the black- hole mass measured by Kollatschny is consistent with our measurement if we use 
FWHM instead of <Tii ne and Kollatschny 's value for the scaling factor /. 

NGC 3227. This AGN was the target of two separate optical campaigns, one by the LAG 
consortium in 1990 (Salamanca et al. 1994) and one at CTIO in 1992 (Winge et al. 1995). The 
rms spectra formed from the LAG data were recently presented by Onken et al. (2003). We 
completely reanalyzed the CTIO data. The original reduced spectra were rescaled in flux using 
the van Groningen & Wanders (1992) algorithm that has been used in most of the International 
AGN Watch campaigns and in the Ohio State program, and new continuum and H/3 emission-line 
light curves were measured from the rescaled spectra. While this resulted in some improvement in 
the H/3 lag determination and uncertainty, the rms spectrum was still quite noisy because of the 
combination of a low amplitude of variability and a relatively insensitive detector (see the notes on 
IC 4392A, below). 

NGC 3516. This is another object observed by the LAG consortium (Wanders et al. 1993). We 
note that the light curves for this object may be less reliable than those of other objects as the 
extended narrow-line region in this object makes narrow line-based flux calibration vulnerable to 
seeing effects, which then have to be modeled (Wanders et al. 1992). We used the scaled and 
corrected spectra to determine the rms spectrum, following Onken et al. (2003). 

NGC 3783. This was the second major multiwavelength campaign undertaken by the Interna- 
tional AGN Watch (Reichert et al. 1994; Stirpe et al. 1994). This object was completely reanalyzed 
by Onken & Peterson (2002); new UV light curves based on IUE NEWSIPS data were measured, 
and the optical spectra were completely recalibrated using the van Groningen & Wanders (1992) 
algorithm. The results presented here are based on the Onken & Peterson reanalysis. 

NGC 4051. This object was also studied recently by Shemmer et al. (2003), who obtain a black 
hole mass consistent with the results of Peterson et al. (2000), who present the data we have used 
here. 

NGC 4151. We analyze two sets of optical data on this object, from the Wise Observatory 
campaign in 1988 (Maoz et al. 1991) and from the International AGN Watch project in 1993-94 
(Kaspi et al. 1996). Unfortunately, the rms line profiles in the 1988 data are too poor to use due to 
a combination of narrow-line residuals in the rms spectra, uncertain narrow-line removal from the 
original spectra, and a variable line-spread function. We therefore do not use these data, except 
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to the extent of noting that the results are broadly consistent with the later AGN Watch results; 
though the lag measurements appear to be reliable, the line widths in the rms spectra are not 
trustworthy. 

PG 1211+143. This AGN is another I Zw 1-like object (see PG 0804+761 above). Our level 
of confidence in the results for this object is low: the lags are suspect because the amplitude of 
variability is low, the variations are slow, and the time sampling is not especially good. The time 
lags are highly uncertain on account of this; the CCFs are very flat-topped and uncertain, as can 
be seen clearly in Fig. 13. Moreover, there are some problems with stability in the Ha region of the 
spectrum that makes the line-width measurements of Ha highly uncertain. Given these difficulties, 
we exclude this object from further analysis. 

PG 1226+023. This is the well-known quasar, 3C 273, which is another I Zw 1-like object 
(see PG 0804+761 above) in which the [Oiii] lines are strongly blended with Fen lines (Peterson, 
Meyers, & Capriotti 1984). The narrow-line components of H7 and H/3 are weak in both the mean 
and rms spectra. However, there are strong narrow-line residuals in the Ha region, so the Ha 
line-width measurements cannot be trusted. 

PG 1229+204. Inspection of the light curves shows that the variations in this object are not 
well sampled. The large differences in lags for the Balmer lines make the results on this object 
rather dubious. We have low confidence in the results for this object. 

NGC 4593. This is another object from the LAG campaign in 1990 (Dietrich et al. 1994), where 
the data have been reanalyzed by Onken et al. (2003). We regard the H/3 lag as completely 
unreliable because it is so much smaller than the mean sampling interval. The Ha lag should also 
be regarded with some caution. 

PG 1307+085. The Ha region of the rms spectrum in this object shows residual narrow-line 
Ha, though this probably affects only the FWHM measurement. 

IC 4329A. This object and NGC 3227 were both observed in the CTIO monitoring program 
(Winge et al. 1995, 1996). These observations employed a Reticon detector, which yielded poorer 
quality spectra than obtained with the CCDs used in virtually every other ground-based campaign. 
We attempted to improve the original lightcurves by rescaling the original spectra in flux by using 
the van Groningen & Wanders (1992) algorithm and remeasuring the continuum and emission-line 
fluxes. This did improve light curves and rms spectra, but only marginally. The light curves are very 
poor and the time-lag measurements should be regarded with caution. The FWHM measurement 
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is very poor (note the large uncertainty yielded by our measurement algorithm). We have little 
confidence in the mass determination for this object. 

Mrk 279. We examined two completely independent sets of data, one from the Wise Observatory 
program in 1988 (Maoz et al. 1990) and one from an International AGN Watch project in 1996 
(Santos-Lleo et al. 2001). The AGN Watch data appear to be quite good. The Wise Observatory 
data, however, have a number of problems: the CCFs for both Ha and H/3 show a correlated error 
signal at zero lag 5 ; this problem is particularly bad at H/3, and we will therefore not use the H/3 
lag measurements. Furthermore, the Ha region of the rms spectrum is strongly contaminated by 
narrow-line residuals, and H/3 is so weak that the FWHM measurement is meaningless. Neverthe- 
less, the mass determination is quite consistent with the results of the AGN Watch program. 

NGC 5548. Some comments on this object appear in section 5. There are more published 
variability data, by far, on this object than any other, and most of the data are very good. Only 
a few problems need to be pointed out. First, the Ha region of the rms spectrum from the Wise 
Observatory campaign (Netzer et al. 1990) has strong narrow-line residuals that render the line- 
width measures unusable. The CCCD for the associated time series is also rather ambiguous, so we 
do not include these measurements in the mass determination. The H/3 profile in the rms spectrum 
from the fifth year of AGN Watch monitoring (i.e., JD2448954-JD2449255) is very unusual, as 
shown in Fig. 14; there are two peaks, a strong central peak and a weak blue peak. The FWHM 
measurement is thus meaningless. Finally, He II A4686 is prominent in the rms spectrum for the 
first year of the AGN Watch campaign (JD2447509-JD2447809), but is too heavily blended with 
the optical Fen blends to measure in the mean spectrum; this line is not included in this part of 
the analysis (cf. Figs. 4 and 5). 

PG 1700+518. The CCCD for this object shows two peaks, one of them clearly ascribable to 
correlated errors in the continuum and emission-line light curves (Fig. 15). Fortunately, the peaks 
are well-separated, and we can exclude the zero-lag correlated error peak from our analysis. This 
significantly increases the H/3 lag and also the black hole mass relative to the original investigation 
by Kaspi et al. (2000). 

3C 390.3. Some of the difficulties with this object have already been discussed in section 5. Part 
of the problem with these data seems to be just the nature of the variability: apart from a large- 
scale outburst during the early part of the monitoring campaign, the variations were very weak. 



5 Correlated errors result from flux calibration problems. An error in flux calibration offsets both the line and 
continuum measurements based on that spectrum in the same direction, thus introducing a spurious cross-correlation 
signal at zero lag. 
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As noted earlier, the lines are very broad and blended, and the most consistent black hole mass is 
obtained by using the line width measured from the unblended side of each line and then assuming 
symmetry. A more detailed attempt at deconvolution might improve this. Finally, we note that in 
this case, we used the UV continuum (at 1370 A) as the driving continuum in the cross-correlation 
analysis. 

Mrk 509. Like Mrk 110, the HeiiA4686 line is prominent in the rms spectrum of this object, 
so we therefore attempted to measure it in each spectrum and produce a light curve. Because of 
blending with optical Fe n emission, the resulting light curve is probably not as reliable as the H/3 
light curve. 

NGC 7469. This object is one of a very few in which a lag has been detected between the UV 
and optical continuum variations (Wanders et al. 1997; Collier et al. 1998). We therefore in this 
case use the UV continuum (at 1315 A) as the driving continuum in the cross-correlation analysis. 

8. THE MASS LUMINOSITY RELATIONSHIP 

Our improved database can be used to investigate the relationships between BLR radius and 
luminosity and black hole mass and luminosity. The radius-luminosity relationship is discussed 
in a companion paper (Kaspi et al. 2004), and here we will discuss only the mass-luminosity 
relationship. 

We computed the optical luminosity from the flux measurements in the original data sources 
that were made over the time intervals used in this analysis. In each case, we selected the continuum 
waveband closest to 5100 A in the AGN rest frame. We corrected for Galactic reddening using the 
extinction values in column (7) of Table 1 and using the reddening curve of Savage & Mathis 
(1979), adjusted to Ay/E(B — V) = 3.2. Luminosity distances were computed using the redshifts 
given in column (5) of Table 1 and by assuming a standard flat ACDM cosmology with S7b = 0.04, 
Qdm = 0.26, S7a = 0.70, and H Q = 70kms _1 Mpc -1 . For AGNs with multiple mass determinations, 
we computed the average value of log \L\ from each individual time series. The optical luminosities 
log \L\ are given in column (4) of Table 8. The corresponding uncertainties represent the amplitude 
of continuum variability during the reverberation experiment. 

In Fig. 16, we plot the reverberation-based masses we have derived as a function of the mean 
luminosity. This figure can be compared directly with Fig. 8 of Kaspi et al. (2000), which reveals 
that a much better-defined mass-luminosity relationship results from our improved analysis. We 
have also estimated the bolometric luminosity in the same fashion as Kaspi et al., i.e., L^ \ ~ 
9AL^(5100 A), and this scale is shown on the top of Fig. 16. The diagonal lines show the Eddington 
limit, and 10% and 1% its value. 
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We note that this constitutes only a preliminary version of the mass-luminosity relationship, 
for comparison with earlier work. A more exhaustive study of this relationship is continuing. 
Several comments are in order: 



1. It is reassuring that there are no objects above the Eddington limit (i.e., to the right of the 
diagonal line), in contrast to the results of Kaspi et al. This is because our analysis has 
corrected a number of errors in earlier work, most notably removal of residual narrow-line H/3 
from many of the rms spectra, which resulted in larger line widths and rather higher overall 
black hole masses. The object closest to the Eddington limit is, not surprisingly, 3C 273 
(PG 1226+023). 

2. The bolometric correction we have assumed is nominal, based on a spectral energy distribution 
(SED) with a strong blue bump. However, observed SEDs suggest that a smaller ratio, e.g., 
-^boi ~ 5AL^ (Netzer 2003), may be on average more appropriate. Moreover, the bolometric 
correction we adopt may not be appropriate for all types of AGNs at arbitrary luminosity. 
The bolometric luminosities are thus uncertain and should be treated with some caution. 

3. The optical luminosities used here have not been corrected for the contribution of starlight 
from the host galaxies. This can be a significant factor, especially in the lower-luminosity ob- 
jects. For example, the standard aperture used in the NGC 5548 optical monitoring program, 
5".0 x 7".5, admits a starlight flux of F x (5100 A) 3.4x10 15 ergs s 1 cm 2 A 1 (Romanishin 
et al. 1995). Correcting for this reduces the optical luminosity entry for NGC 5548 in Table 
8 by ~ 0.24 dex. A program is currently underway to determine the starlight contribution to 
the optical luminosity of each of these objects. 

4. Internal extinction has not been accounted for in any way. Correction for extinction will 
move objects to the right in this diagram. It is worth noting that the object with the lowest 
Eddington rate (farthest to the left of the O.OlLEdd diagonal) is NGC 3227, which is a rather 
dusty object (Pogge & Martini 2002). The low luminosity of this object relative to its black 
hole mass may be a result of internal extinction. 

5. Given the small formal error bars for most of the objects, we believe that much of the scatter 
in Fig. 16 is real. Moreoever, we find that the scatter correlates with other AGN properties; 
the I Zw 1-type objects generally lie along the bottom edge of the envelope defined by the data 
points. In Fig. 16, the NLSls are shown as open circles, and all of them except NGC 4051 lie 
on the lower edge of the mass-luminosity envelope. Conversely, the one object with strongly 
double-peaked Balmer line profiles, 3C 390.3, lies along the upper edge of the envelope. The 
locations of these extreme objects on this diagram suggest that at least some of the dispersion 
of the data points correlates with Eigenvector 1, consistent with the suggestion originated 
by Boroson h Green (1992) and reaffirmed by numerous later authors that Eigenvector 1 
appears to be driven by Eddington ratio rh = M/MEdd- However, the physical origin of the 
scatter observed in Fig. 16 could be attributable either to differences in Eddington ratio or to 
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inclination effects. Decreasing inclination (i.e., from edge-on to face-on) will translate points 
to the right as the apparent luminosity increases on account of decreased limb darkening and 
downward as the rotational velocities appear to decrease. Increasing the Eddington ratio will 
translate points in the same sense. 

6. The best fit mass-luminosity relationship is found to be 

log(M/10 8 M Q ) = -0.12(±0.07) +0.79(±0.09)log (AL A (5100 A)/10 44 ergs" 1 ) . (9) 

However, there is no reason to believe that there are no selection effects operating. Interest- 
ingly, the lower edge of the envelope seems to parallel the lines of constant Eddingtion ratio 
rather well, suggesting that the intrinsic mass-luminosity slope may not differ signficantly 
from unity. 

9. SUMMARY 

In this contribution, we have improved the calibration of the reverberation-based AGN black 
hole mass scale by decreasing random and systematic errors and by drawing on the AGN Mbh _ c* 
relationship to establish a statistically accurate calibration that is tied to other methods of black 
hole mass measurement. We have undertaken a consistent reanalysis of a highly inhomogeneous 
database that consists of 117 separate time series, not including several others that were deemed to 
be too poor to use, in the process accounting for a variety of systematic effects such as time dilation 
on the time lags and spectral resolution on the line widths. Each time series is treated independently 
and yields an independent estimate of the black hole mass, thus reducing random uncertainties. 
Poor or suspicious data are removed from the database, noting especially the susceptibility of time- 
lag measurements to correlated errors and other types of aliasing. We find that the most consistent 
mass measurements are obtained by using the cross-correlation centroid r ccnt to characterize the 
light-travel time across the BLR and by using the line dispersion an ne as measured in the rms 
spectrum to characterize the velocity dispersion of the BLR. In practice, special care has been taken 
to remove residual narrow-line contamination of the rms spectrum in cases where it is present. The 
result of this analysis is a revised AGN mass scale based on 35 reverberation-mapped AGNs that 
is statistically accurate to better than a factor of three. 
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A. SOME COMMENTS ON CROSS-CORRELATION METHODOLOGY 

Welsh (1999) has suggested a number of modifications to our error analysis procedures and we 
have investigated these through additional Monte Carlo simulations similar to those described by 
P98b. 

Emission-line lag measurements are made by cross-correlating the emission-line and continuum 
light curves. Since the data are almost never regularly sampled, real data points in one time series 
are matched with values obtained by linear interpolation of the other time series. The cross- 
correlation function (CCF) is measured twice, once interpolating in the continuum series, and once 
in the emission-line series, and the final CCF is determined by averaging these two results. The 
CCF is characterized by (1) its peak value (highest value of the correlation coefficient) r max , (2) 
the time delay corresponding to this value r pca k, and (3) the centroid r ccnt of the peak in the CCF. 
In practice, r cent is evaluated using only points with values about some threshold, usually 0.8r max . 

To assess the uncertainties in the determination of Tp ea k and T cent , we use the model-independent 
Monte Carlo FR/RSS method described by P98b. For each Monte Carlo realization, we start with 
a parent light curve of N data points and from this make N independent random selections of these 
data points without consideration for previous selections. The redundant selections are then dis- 
carded, leaving a new light curve of M < N points; typically, the fraction of points from the parent 
light curve that remain unselected in each realization is ~ 1/e. We refer to this process as "random 
subset selection" (RSS), and it seems to successfully account for uncertainties due to effects of 
individual data points. Another major source of uncertainty is the uncertainty in the measured 
continuum and emission-line fluxes, and these can be significant if the amplitude of variability is 
not much larger than the flux errors in individual data points. We attempt to account for this 
by altering the fluxes of the M data points in each Monte Carlo realization by random Gaussian 
deviates scaled by the flux error associated with the data point. We refer to this process as "flux 
randomization" (FR). This process is repeated for many independent realizations (usually 2000 or 
more for the light curves analyzed here), and for the CCF from each realization, r max , T pca k, and 
r C ent are recorded. These are used to build up "cross-correlation peak distributions (CCPDs)" and 
"cross-correlation centroid distributions (CCCDs)" for r pca k and T cen t> respectively. This is done 
because the distributions of these values are rarely even approximately Gaussian (Maoz & Netzer 
1989). 

Welsh (1999) points out a number of potential problems with cross-correlation methodology, 
some of which can produce biases in determination of the emission-line lags. However, simulations 
that mimic as much as possible real AGN behavior have not revealed any strong systematic biases. 
Nevertheless, we caution that cross-correlation of light curves to measure emission- line lags is a 
rather crude tool, but one that seems to be effective with the limited AGN variability data available 
at present. Welsh has also suggested a number of modification to the FR/RSS method, and we 
have extended the Monte Carlo simulations described by P98b to test these suggestions under what 
we regard as reasonably realistic conditions, and we describe our results below. 
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Flux Uncertainties and Redundant Selections. Welsh (1999) suggests an alternative strat- 
egy to FR/RSS, namely counting redundant selections in each realization and appropriately reduc- 
ing the flux error for multiply-selected points. Specifically, the flux uncertainty for each point that 
is selected 1 < n < N times should be reduced by a factor of n 1 / 2 . Thus, rather than omitting 
redundant points, each of the N selections has some real effect on the outcome of the realization and 
the weighting is more consistent with the standard bootstrap method (e.g., Diaconis & Efron 1983) 
on which the FR/RSS method is based. We carried out detailed simulations like those described 
by P98b to determine the efficacy of Welsh's proposed scheme. We compare different error assess- 
ments by examining the width of the CCPDs and CCCDs produced by simulations; we presume 
that the algorithm that yields the narrowest CCPDs and CCCDs (i.e., the highest precision mea- 
surements) is the best, as long as the errors are not underestimated. Our simulations confirm that 
the error estimates using Welsh's method are superior to those of the original FR/RSS algorithm. 
The uncertainties in r pea k are typically lower by ~ 8%, and the uncertainties in r ccn t are lower by 
~ 3%. Following the procedures of P98b, we have also carried out model-dependent Monte Carlo 
simulations in which we use a known model for the transfer function in order to verify that Welsh's 
algorithm does not underestimate the uncertainties. We therefore adopt this improvement in the 
cross-correlation analysis employed in this contribution. 

Detrending the Light Curves. Welsh (1999) also suggests that cross-correlation results are 
more accurate if the data are first "detrended," i.e., a low-order (usually linear) polynomial is 
fitted to the light curve and subtracted off prior to carrying out the cross-correlation analysis. Our 
simulations support this, but only in the case where the sampling is excellent, i.e., long duration 
at high resolution (for the cases considered by P98b, say, a 200-day experiment with observations 
once per day). However, we find that under conditions of more marginal sampling (say, a 200-day 
experiment with only 40 observations), which includes nearly all reverberation data that exist, 
detrending either makes little difference or can, in fact, lead to occasional gross errors in the lag 
determination and, consequently, gross overestimates of the uncertainties in the lags. For this 
reason, we elect to not detrend our data prior to cross-correlation. 

Peak or Centroid? Whether or not the cross-correlation lag is better characterized by the peak 
or centroid of the CCF has been debated on many occasions. The advantage of the centroid is that 
it is related to the centroid of the transfer function and is better defined when the CCF has a broad 
peak, and for these reasons, we prefer it. Welsh's simulations, however, suggest that the peak is a 
more robust measure than the centroid, quite at odds with what was found by P98b and contrary 
to our general experience. The difference seems to be attributable to differences in the transfer 
functions used by Welsh and by P98b; P98b used only thin-shell and thick-shell transfer functions 
which are flat-topped (i.e., with poorly defined peaks), whereas Welsh used only Gaussian transfer 
functions. If we repeat the simulations of P98b with Gaussian transfer functions, we find little 
reason to prefer one measure to the other. Further investigation of this issue using real data, as 
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described in this paper, leads us to continue to prefer the CCF centroid, although the CCD peak 
is also an acceptable way to characterize the emission-line lags. 

Centroid Threshold. Aliasing effects can lead to complex structures in CCFs; rarely is the 
principal peak in the CCF isolated and well-defined. It is therefore necessary to compute the CCF 
centroid only over a restricted range, including points with values larger than some fixed fraction of 
the peak value r max . It is conventional to use a threshold of 0.8r max for this computation, although 
in some cases where the peak is noisy, lower thresholds are used. In some cases, the centroid can 
vary significantly for different threshold selections (Koratkar & Gaskell 1991), though this appears 
to be less of a problem with well-sampled data (e.g., Dietrich et al. 1993). Simulations based on 
those of P98b, however, show that a threshold value of 0.8r max is generally a good choice. With 
lower thresholds, we find that the CCCD is significantly broadened. For a threshold of 0.5r max , 
for example, we find that the width of the CCCD increases by 10-20%, depending somewhat on 
the details of the transfer function (i.e., more sharply peaked transfer functions give results less 
sensitive to the selected threshold value, as one might expect). We will therefore continue to use a 
threshold of 0.8r max , unless otherwise noted. 
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Fig. 1. — Measurement of FWHM for double-peaked profiles. A continuum is interpolated un- 
derneath the line profile and is subtracted from the data. Peak fluxes are identified on the 
short-wavelength and long-wavelength peaks, P(X)max,short and P(X)max,iong, respectively. On 
each side of the line, wavelengths corresponding to the half-maximum fluxes 0-5P(X) rnax ,short and 
0.5P(A) 

max,long are found moving upward from the continuum at Ai and downward from the peak 
at A2. In this example, Ai = A2 on the short-wavelength side only. The nominal position of the 
half- maximum point is taken to be the average of Ai and A2. 
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Fig. 2. — Comparison of FWHM measurements and uncertainties from this paper with those from 
WPM. For the sake of meaningful comparison, the new measurements here have not been corrected 
for spectral resolution, as have all other line-width measurements in this paper. 



-32- 




Fig. 3. — Line widths versus time lags for emission lines in the rms spectra of NGC 5548. The top 
row shows the line dispersion <7ii ne as a line width measure, and the bottom row shows FWHM. 
The left column shows the CCF centroid r cen t as the time-lag measure and the right column shows 
the location of the CCF peak r pea k. The solid line is the best fit to the data, and the dotted line 
is a forced fit to slope —1/2, the virial slope. The fit parameters are summarized in Table 7. The 
open circles are measurements of H/3 for 14 different years. The filled circles represent all of the 
other lines. 
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Fig. 4. — Line width versus time lags for emission lines in the mean spectra of NGC 5548. The 
data are plotted as in Fig. 3. Note the dramatic change in the scatter, especially for the H/3 line, 
by using the mean rather than the rms line profile for the line-width measurement. 
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Fig. 5. — Line widths versus time lags for emission lines in the mean spectra of NGC 5548, but 
with the narrow component of H/3 removed. The data are plotted as in Figs. 3 and 4. Removing 
the narrow-component of H/3 greatly reduces the scatter around the virial relationship, but the 
scatter is still much more pronounced than in Fig. 3, where the line- widths are measured in the rms 
profiles. The measurement uncertainties on o"ii ne are too small to show up on this diagram. Note 
that the fits to the data in the lower panel are unweighted fits on account of the large number of 
outliers. 
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Fig. 6. — Line widths versus time lags for emission lines in the rms spectra of NGC 3783. The data 
are plotted as in Fig. 3. 
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Fig. 7. — Line widths versus time lags for emission lines in the rms spectra of NGC 7469. The data 
are plotted as in Fig. 3. 
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Fig. 8. — Line width versus CCF centroid r ceil t for emission lines in the rms spectra of 3C 390.3. 
The top row shows the line dispersion <Tii ne used as the line-width measure, and the bottom row 
shows FWHM. The left-hand column shows our standard line-width measurement, using the whole 
line, but with no effort to account for blending. In the right-hand column, we attempt to account 
for blending by measuring the unblended half of each line. In any case, the large uncertainties in 
the time lags preclude a critical test of a virial relationship. 
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Fig. 9. — Distribution of values of FWHM/cij ne for all lines in the rms spectra used in this analysis 
(i.e., columns (5) and (6) of Table 6), excluding highly uncertain values (preceeded with a colon in 
Table 6). The mean and standard deviation for the total distribution are 2.03 ± 0.59. The black 
area shows the distribution for multiple measurements of H/3 in NGC 5548 (the mean and standard 
deviation for this subset are 2.73 ± 0.24). The gray area shows the distribution for other lines in 
NGC 5548 (mean and standard deviation for all lines in NGC 5548, including H/3, are 2.45 db 0.44. 
The vertical dotted line is at FWHM/o"ii ne = 2.355, which is appropriate for a Gaussian line 
profile. Values smaller than this indicate lines that have weaker cores and strong wings relative to 
a Gaussian. 
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Fig. 10. — Distribution of values of T pca k/T ccnt for all lines in this analysis (i.e., columns (3) and 
(4) of Table 6), excluding highly uncertain values (preceeded with a colon in Table 6). The mean 
and standard deviation for this distribution are 0.95 ± 0.20. 
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Fig. 11. — Cross-correlation centroid distribution for the continuum-H/3 cross-correlation for Mrk 
79 during the period JD2449996 to JD2450220. It is not obvious which peak corresponds to the 
correct lag. 
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Fig. 12. — Mean (top panel) and rms (bottom panel) spectra of Mrk 110. The rms spectrum shows 
clearly the very broad He n A4686 and broad H/3. The narrow [O III] AA4959, 5007 lines appear only 
in the mean spectrum, except for weak residuals in [O m] A5007 that appear in the rms spectrum. 



-42 - 




Fig. 13. — Cross-correlation functions for PG 1211+143. The dotted line is the Ha CCF, the solid 
line is the H/3 CCF, and the dashed line is the H7 CCF. The centroids and peaks are poorly defined 
because the CCFs are broad and flat-topped. 
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Fig. 14. — Mean (top panel) and rms (bottom panel) H/3-region spectra of NGC 5548 during the 
period JD2448954 to JD2449255, the fifth year of monitoring by the International AGN Watch 
(1993). The H/3 profile in the rms spectrum is double-peaked; there is a strong peak at line center 
and another in the shortward wing of the line. 
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Fig. 15. — Cross-correlation centroid distribution for the continuum-H/3 cross-correlation for PG 
1700+518. The peak at zero lag is clearly ascribable to correlated error, so it can be rejected. 
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Fig. 16. — Black hole mass vs. luminosity for 35 reverberation mapped AGNs. The luminosity 
scale on the lower x-axis is logAL^ in units of ergs s _1 . The upper x-axis shows the bolometric 
luminosity assuming that L^ \ k, 9\L\. The diagonal lines show the Eddington limit L^d, 0.1I/Edd> 
and O.OlLEdd- The open circles represent NLSls. Other labeled points are discussed in the text. 
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12 


10 


.49 


0.02578 


13.7 


0.153 


J000619. 5+201211 


PG 0003+199 


PG 0026+129 


2 


00 


29 


13.600 


+13 


16 


03 


.00 


0.14200 


15.4 


0.307 


J002913. 8+131605 




PG 0052+251 


2 


00 


54 


52.100 


+25 


25 


38 


.00 


0.15500 


15.4 


0.205 


J005452. 2+252539 




Fairall 9 


3,4 


01 


23 


45.780 


—58 


48 


20 


.50 


0.04702 


13.5 


0.116 


J012345. 8-584821 




Mrk 590 


1 


02 


14 


33.562 


—00 


46 


00 


.09 


0.02638 


13.8 


0.161 


J021433. 6-004600 


NGC 863 


3C 120 


1 


04 


33 


11.095 


+05 


21 


15 


.62 


0.03301 


14.1 


1.283 


J043311. 1+052115 


Mrk 1506 


Akn 120 


1 


05 


16 


11.421 


—00 


08 


59 


.38 


0.03230 


14.1 


0.554 


J051611. 4-000900 


Mrk 1095 


Mrk 79 


1 


07 


42 


32.797 


+49 


48 


34. 


.75 


0.02219 


13.9 


0.305 


J074232. 8+494835 




PG 0804+761 


2 


08 


10 


58.600 


+76 


02 


42 


.00 


0.10000 


15.1 


0.150 


J081058. 5+760243 




PG 0844+349 


2 


08 


47 


42.400 


+34 


45 


04 


.00 


0.06400 


14.0 


0.159 


J084742.5+344505 




Mrk 110 


1 


09 


25 


12.870 


+52 


17 


10 


.52 


0.03529 


16.0 


0.056 


J092512. 9+521711 




PG 0953+414 


2 


09 


56 


52.400 


+41 


15 


22 


.00 


0.23410 


14.5 


0.054 


J095652. 3+411522 




NGC 3227 


5,6,7 


10 


23 


30.589 


+19 


51 


53 


.99 


0.00386 


11.1 


0.098 


J102330. 6+195156 




NGC 3516 


6,8 


11 


06 


47.490 


+72 


34 


06 


.88 


0.00884 


12.5 


0.183 


J110647. 4+723407 




NGC 3783 


9,10,11 


11 


39 


01.77 


—37 


44 


18 


.7 


0.00973 


12.6 


0.514 


J113901. 8-374419 




NGC 4051 


12 


12 


03 


09.614 


+44 


31 


52 


.80 


0.00234 


10.8 


0.056 


J120309. 6+443153 




NGC 4151 


13,14 


12 


10 


32.579 


+39 


24 


20 


.63 


0.00332 


11.5 


0.119 


J121032. 5+392421 




PG 1211+143 


2 


12 


14 


17.700 


+14 


03 


12 


.60 


0.08090 


14.6 


0.150 


J121417. 7+140313 




PG 1226+023 


2 


12 


29 


06.700 


+02 


03 


08 


.60 


0.15834 


12.8 


0.089 


J122906. 7+020308 


3C 273 


PG 1229+204 


2 


12 


32 


03.605 


+20 


09 


29 


.21 


0.06301 


15.3 


0.117 


J123203. 6+200930 


Mrk 771 
Ton 1542 


NGC 4593 


6,15 


12 


39 


39.425 


-05 


20 


39 


.34 


0.00900 


11.7 


0.106 


J123939. 4-052039 


Mrk 1330 


PG 1307+085 


2 


13 


09 


47.000 


+08 


19 


48. 


.90 


0.15500 


15.3 


0.145 


J130947.0+081949 




IC 4329A 


16 


13 


49 


19.29 


-30 


18 


34. 


.4 


0.01605 


14.0 


0.255 


J134919. 3-301834 




Mrk 279 


17,18 


13 


53 


03.447 


+69 


18 


29 


.57 


0.03045 


14.6 


0.068 


J135303. 5+691830 




PG 1351+640 


2 


13 


53 


15.808 


+63 


45 


45 


.41 


0.08820 


14.8 


0.088 


J135315. 7+634546 




PG 1411+442 


2 


14 


13 


48.300 


+44 


00 


14 


.00 


0.08960 


15.0 


0.036 


J141348. 3+440014 




NGC 5548 


19,20,21,22,23 


14 


17 


59.534 


+25 


08 


12 


.44 


0.01717 


13.3 


0.088 


J141759. 6+250813 




PG 1426+015 


2 


14 


29 


06.588 


+01 


17 


06. 


.48 


0.08647 


17.5 


0.137 


J142906. 5+011704 


Mrk 1383 


Mrk 817 


1 


14 


36 


22.068 


+58 


47 


39 


.38 


0.03145 


14.5 


0.029 


J143622. 1+584740 


PG 1434+590 


PG 1613+658 


2 


16 


13 


57.179 


+65 


43 


09 


.58 


0.12900 


15.2 


0.114 


J161357.2+654309 


Mrk 876 


PG 1617+175 


2 


16 


20 


11.288 


+17 


24 


27 


.70 


0.11244 


15.5 


0.180 


J162011. 3+172428 


Mrk 877 


PG 1700+518 


2 


17 


01 


24.800 


+51 


49 


20 


.00 


0.29200 


15.4 


0.151 


J170125. 0+514920 




PG 1704+608 


2 


17 


04 


41.370 


+60 


44 


30 


.50 


0.37100 


15.3 


0.097 


J170441. 5+604428 




3C390.3 


24,25 


18 


42 


08.9899 


+79 


46 


17 


.127 


0.05610 


14.4 


0.308 


J184209. 0+794617 




Mrk 509 


1 


20 


44 


09.738 


-10 


43 


24. 


.54 


0.03440 


13.0 


0.248 


J204409. 7-104324 




PG 2130+099 


2 


21 


32 


27.813 


+ 10 


08 


19 


.46 


0.06298 


14.6 


0.192 


J213227.8+100819 


II Zw 136 

Mrk 1513 


NGC 7469 


26,27 


23 


03 


15.623 


+08 


52 


26 


.39 


0.01632 


13.0 


0.297 


J230315.6+085226 


Mrk 1514 



References: 1: Peterson et al. 1998a; 2: Kaspi et al. 2000; 3: Santos-Lleo et al. 1997; 4: Rodrigucz-Pascual ct al. 1997; 5: Salamanca et al. 
1994; 6: Onken ct al. 2003; 7: Winge et al. 1995; 8: Wanders et al. 1993; 9: Stirpc et al. 1994; 10: Onken & Peterson 2002; 11: Reichert et 
al. 1994; 12: Peterson ct al. 2000; 13: Kaspi ct al. 1996; 14: Maoz et al. 1991; 15: Dietrich ct al. 1994; 16: Winge et al. 1996; 17: Santos-Lleo 
ct al. 2001; 18: Maoz et al. 1990; 19: Peterson et al. 2002 and references therein; 20: Dietrich et al. 1993; 21: Clavcl ct al. 1991; 22: Korista 
ct al. 1995; 23: Netzer ct al. 1990; 24: Dietrich et al. 1998; 25: O'Brien et al. 1998; 26: Collier et al. 1998; 27: Wanders et al. 1997. 
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Table 2. Wavelength Windows for Line Width Measurements 



Data Res. Julian Dates Continuum Continuum Line 
Object Ref. 1 Set (A) (-2400000) Window (A) Window (A) Limits (A) 

(1) (2) (3) (4) (5) (6) (7) (8) 



Mrk 335 


1 


H/3 


9.6 


49156-49338 


4900-4910 


5095-5105 


4930-5030 




1 


H/3 


9.6 


49889-50118 


4910-4920 


5100-5110 


4930-5030 


PC 0026+129 


2 


Ha 


12.5 


48836-51084 


7290-7300 


7700-7740 


7304-7728 




2 


H/3 


12.5 


48545-51084 


5442-5450 


5820-5840 


5460-5700 


PC 0052+251 


2 


Ha 


12.5 


48837-51084 


7028-7108 


7725-7800 


7398-7724 




2 


H/3 


12.5 


48461-51084 


5418-5488 


5908-5978 


5518-5698 




2 


H 7 


12.5 


48461-51084 


4864-4914 


5418-5488 


4926-5118 


Fairall 9 


3 


H/3 


11.2 


50473-50665 


4936-4946 


5295-5323 


4946-5295 




4 


Civ A1549 


6.0 


50473-50713 


1540-1565 


1760-1770 


1565-1670 




4 


Lya 


6.0 


50473-50713 


1227-1232 


1370-1380 


1232-1334 


Mrk 590 


1 


H/3 


9.6 


48090-48323 


4940-4955 


5050-5065 


4965-5045 




1 


H/3 


9.6 


48848-49048 


4880-4890 


5235-5245 


4890-5095 




1 


H/3 


9.6 


49183-49338 


4930-4935 


5165-5170 


4955-5069 




1 


H/3 


9.6 


49958-50122 


4920-4930 


5170-5180 


4930-5068 


3C 120 


1 


H/3 


9.6 


47837-50388 


4958-4972 


5100-5108 


4972-5100 


Akn 120 


1 


H/3 


9.6 


48148-48344 


4910-4920 


5246-5255 


4920-5107 




1 


H/3 


9.6 


49980-50175 


4910-4920 


5245-5255 


4920-5102 


Mrk 79 


1 


H/3 


9.6 


47838-48044 


4795-4810 


5190-5210 


4850-5040 




1 


H/3 


9.6 


48193-48393 


4850-4870 


5075-5085 


4890-5040 




1 


H/3 


9.6 


48905-49135 


4880-4885 


5080-5086 


4885-5040 




1 


H/3 


9.6 


49996-50220 


4830-4840 


5050-5060 


4840-5040 


PC 0804+761 


2 


Ha 


12.5 


48319-51085 


6746-6814 


7472-7520 


6940-7412 




2 


H/3 


12.5 


48319-51085 


5224-5264 


5598-5630 


5266-5474 




2 


H7 


12.5 


48319-51085 


4680-4705 


4880-4905 


4710-4850 


PC 0844+349 


2 


Ha 


12.5 


48319-51085 


6708-6764 


7248-7348 


6848-7120 




2 


H/3 


12.5 


48319-51085 


5048-5078 


5414-5440 


5098-5248 




2 


H7 


12.5 


48319-51085 


4503-4535 


4710-4750 


4536-4690 


Mrk 110 


1 


H/3 


9.6 


48953-49149 


4970-4980 


5090-5100 


4980-5090 




1 


H/3 


9.6 


49751-49874 


4715-4740 


5250-5260 


4975-5066 




1 


H/3 


9.6 


50010-50262 


4965-4970 


5150-5160 


4975-5070 




1 


HeiiA4686 


9.6 


48953-49149 


4730-4740 


4970-4980 


4740-4970 




1 


HeiiA4686 


9.6 


49751-49874 


4740-4750 


5240-5260 


4750-4965 




1 


HeiiA4686 


9.6 


50010-50262 


4750-4760 


4960-4970 


4760-4960 


PC 0953+414 


2 


H/3 


12.5 


48319-50997 


5846-5892 


6288-6320 


5909-6086 




2 


H 7 


12.5 


48319-50997 


5238-5278 


5472-5506 


5284-5450 


NGC 3227 


5,6 


Ha 


9.1 


47894-48045 


6420-6430 


6800-6830 


6430-6780 




5,6 


H/3 


9.1 


47894-48049 


4780-4795 


5115-5120 


4800-4955 




7 


H/3 


11.2 


48623-48776 


4765-4785 


5085-5100 


4800-4940 



-48- 



Table 2 — Continued 



Data Res. Julian Dates Continuum Continuum Line 

Object Ref. 1 Set (A) (-2400000) Window (A) Window (A) Limits (A) 

(1) (2) (3) (4) (5) (6) (7) (8) 



NGC 3516 


6,8 


Ha 


9.1 


47894-48047 


6426-6477 


6820- 


-6850 


6477- 


-6820 




6,8 


H/3 


9.1 


47894-48047 


4782-4807 


5175- 


-5195 


4807- 


-5031 


NGC 3783 


9,10 


H/3 


11.2 


48607-48833 


4800-4810 


5128- 


-5160 


4810- 


-4984 




10,11 


SiivA1400 


6.0 


48611-48833 


1363-1373 


1464- 


-1470 


1373- 


-1450 




10,11 


CivA1549 


6.0 


48611-48833 


1464-1470 


1710- 


-1725 


1509- 


-1613 




10,11 


HenA1640 


6.0 


48611-48833 


1464-1470 


1710- 


-1725 


1613- 


-1709 


NGC 4051 


12 


H/3 


9.6 


50183-50263 


4810-4820 


4915- 


-4930 


4835- 


-4910 




12 


HenA4686 


9.6 


50183-50263 


4570-4600 


4745- 


-4760 


4610- 


-4740 


NGC 4151 


13 


H/3 


12.5 


49305-49404 


4785-4795 


4946- 


-4954 


4795- 


-4945 




13 


Ha 


12.5 


49305-49404 


6200-6250 


6790- 


-6810 


6360- 


-6737 




14 


H/3 


13.3 


47145-47360 


4767-4775 


5090- 


-5100 


4775- 


-4923 




14 


Ha 


13.3 


47145-47360 


6400-6430 


6682- 


-6690 


6440- 


-6682 


PC 1211+143 


2 


Ha 


12.5 


48319-51000 


6670-6720 


7480- 


-7540 


6900- 


-7310 




2 


H/3 


12.5 


48319-51000 


5120-5160 


5490- 


-5530 


5180- 


-5340 




2 


H 7 


12.5 


48319-51000 


4525-4560 


4780- 


-4810 


4625- 


-4755 


PG 1226+023 


2 


Ha 


12.5 


48361-50997 


7034-7085 


7789- 


-7808 


7250- 


-7788 




2 


H/3 


12.5 


48361-50997 


5494-5522 


5888- 


-5928 


5530- 


-5756 




2 


H 7 


12.5 


48361-50997 


4925-4945 


5135- 


-5165 


4946- 


-5130 


PG 1229+204 


2 


Ha 


12.5 


48319-50997 


6562-6654 


7244- 


-7320 


6838- 


-7148 




2 


H/3 


12.5 


48319-50997 


5052-5088 


5412- 


-5456 


5094- 


-5244 




2 


H 7 


12.5 


48986-50997 


4500-4560 


4694- 


-4758 


4560- 


-4690 


NGC 4593 


6,15 


Ha 


9.1 


47894-48049 


6450-6470 


6735- 


-6760 


6510- 


-6710 




6,15 


H/3 


9.1 


47894-48049 


4745-4790 


5105- 


-5130 


4828- 


-4961 


PG 1307+085 


2 


Ha 


12.5 


49130-51000 


7050-7100 


7750- 


-7800 


7400- 


-7740 




2 


H/3 


12.5 


48319-51042 


5458-5488 


5864- 


-5914 


5508- 


-5698 




2 


H 7 


12.5 


48319-51042 


4870-4915 


5104- 


-5154 


4938- 


-5088 


IC 4329A 


16 


H/3 


11.5 


48643-48832 


4780-4850 


5160- 


-5200 


4850- 


-5060 


Mrk 279 


17 


H/3 


9.6 


50095-50289 


4940-4950 


5185- 


-5200 


4950- 


-5090 




18 


H/3 


13.3 


47205-47360 


4840-4880 


5270- 


-5300 


4950- 


-5090 




18 


Ha 


13.3 


47205-47360 


6520-6600 


6878- 


-6882 


6688- 


-6878 


PG 1411+442 


2 


Ha 


12.5 


48319-51038 


6700-6760 


7490- 


-7560 


6870- 


-7360 




2 


H/3 


12.5 


48319-51038 


5155-5185 


5520- 


-5555 


5201- 


-5382 


NGC 5548 


19 


H/3 


9.6 


47509-47809 


4840-4860 


5030- 


-5060 


4860- 


-5026 




19 


H/3 


9.6 


47861-48179 


4820-4840 


5030- 


-5060 


4840- 


-5030 




19 


H/3 


9.6 


48225-48534 


4840-4850 


5030- 


-5035 


4850- 


-5025 




19 


H/3 


9.6 


48623-48898 


4840-4850 


5030- 


-5035 


4850- 


-5025 




19 


H/3 


9.6 


48954-49255 


4820-4835 


5020- 


-5035 


4840- 


-5000 




19 


H/3 


9.6 


49309-49636 


4760-4770 


5116- 


-5140 


4783- 


-5026 



Table 2 — Continued 



Data Res. Julian Dates Continuum Continuum Line 

Object Ref. 1 Set (A) (-2400000) Window (A) Window (A) Limits (A) 

(1) (2) (3) (4) (5) (6) (7) (8) 





19 


H/3 


9.6 


49679- 


-50008 


4810-4830 


5040-5060 


4830-5010 




19 


H/3 


9.6 


50044- 


-50373 


4845-4855 


5060-5080 


4860-5030 




19 


H/3 


9.6 


50434- 


-50729 


4840-4850 


5025-5044 


4850-5020 




19 


H/3 


9.6 


50775- 


-51085 


4850-4865 


5050-5070 


4865-5030 




19 


H/3 


9.6 


51142- 


-51456 


4860-4875 


5050-5060 


4875-5030 




19 


H/3 


9.6 


51517- 


-51791 


4830-4850 


5050-5060 


4850-5050 




19 


H/3 


9.6 


51878- 


-52174 


4825-4840 


5050-5060 


4840-5050 




20 


HeiiA4686 


9.6 


47509- 


-47797 


4600-4620 


4850-4862 


4620-4850 




21 


SiivA1400 


6.0 


47510- 


-47745 


1380-1390 


1455-1460 


1390-1455 




21 


CivA1549 


6.0 


47510- 


-47745 


1490-1500 


1730-1735 


1500-1630 




21 


HeiiA1640 


6.0 


47510- 


-47745 


1490-1500 


1730-1735 


1630-1730 




21 


Cm] A1909 


6.0 


47510- 


-47745 


1890-1895 


1970-1975 


1895-1970 




22 


SiivA1400 


1.9 


49097- 


-49135 


1370-1380 


1480-1490 


1380-1480 




22 


CivA1549 


1.9 


49060- 


-49135 


1470-1480 


1730-1735 


1525-1615 




22 


HeiiA1640 


1.9 


49097-49135 


1470-1480 


1730-1735 


1615-1730 




22 


Cm] A1909 


1.9 


49097- 


-49135 


1860-1875 


1990-2000 


1875-1990 




23 


H/3 


13.3 


47212- 


-47360 


4820-4825 


5186-5200 


4840-5040 




23 


Ha 


13.3 


47212- 


-47360 


6500-6550 


6770-6780 


6570-6770 


PC 1426+015 


2 


Ha 


12.5 


48334- 


-51042 


6680-6730 


7490-7550 


6800-7450 




2 


H/3 


12.5 


48334- 


-51042 


5105-5145 


5515-5555 


5160-5500 


Mrk 817 


1 


H/3 


9.6 


49000- 


-49212 


4905-4915 


5240-5270 


4915-5073 




1 


H/3 


9.6 


49404- 


-49528 


4905-4925 


5120-5130 


4930-5110 




1 


H/3 


9.6 


49752- 


-49924 


4905-4925 


5120-5130 


4920-5075 


PC 1613+658 


2 


H/3 


12.5 


48397- 


-51073 


5300-5350 


5740-5790 


5360-5566 


PC 1617+175 


2 


Ha 


12.5 


48362- 


-51085 


6898-6968 


7540-7580 


7118-7498 




2 


H/3 


12.5 


48362- 


-51085 


5250-5295 


5672-5722 


5310-5550 


PC 1700+518 


2 


H/3 


12.5 


48378- 


-51084 


6060-6080 


6400-6440 


6200-6400 


3C 390.3 


24 


H/3 


9.6 


49718- 


-50012 


5000-5006 


5246-5255 


5006-5246 




24 


HeiiA4686 


9.6 


49718- 


-50012 


4820-4825 


5000-5006 


4825-5000 




25 


Lya+Nv 


6.0 


49718- 


-50147 


1230-1236 


1345-1350 


1236-1345 




25 


CivA1549 


6.0 


49718- 


-50147 


1522-1535 


1740-1760 


1570-1716 


Mrk 509 


1 


H/3 


9.6 


47653- 


-50374 


4945-4950 


5205-5220 


4950-5087 




1 


HeiiA4686 


9.6 


47653- 


-50374 


4693-4710 


4945-4950 


4710-4945 


PC 2130+099 


2 


Ha 


12.5 


48382- 


-51084 


6760-6790 


7200-7230 


6800-7200 




2 


H/3 


12.5 


48382- 


-51084 


5070-5080 


5440-5450 


5080-5274 




2 


H7 


12.5 


48382- 


-51084 


4510-4540 


4716-4735 


4540-4716 


NGC 7469 


26 


H/3 


12.5 


50237- 


-50295 


4855-4860 


5180-5200 


4870-5014 




26 


Ha 


12.5 


50237- 


-50295 


6500-6530 


6900-6950 


6550-6800 



-50- 



Table 2 — Continued 



Data Res. Julian Dates Continuum Continuum Line 
Object Ref. 1 Set (A) (-2400000) Window (A) Window (A) Limits (A) 

(1) (2) (3) (4) (5) (6) (7) (8) 



27 SiivA1400 6 50245-50293 
27 CivA1549 6 50245-50293 
27 HeiiA1640 6 50245-50293 



1360-1375 1466-1485 1375-1466 
1526-1538 1700-1715 1538-1617 
1526-1538 1700-1715 1617-1700 



1 References as in Table 1. 



Table 3. [O in] A5007 Line Widths 





FWHMQOm] A5007) 


Object 


(km s" 1 ) 


(1) 


(2) 


Mrk 335 


280 


Fairall 9 


425 


Mrk 590 


400 


Akn 120 


490 


Mrk 79 


350 


NGC 3227 


485 


NGC 3516 


250 


NGC 3783 


230 


NGC 4051 


190 


NGC 4151 


425 


NGC 4593 


255 


IC 4329A 


550 


Mrk 279 


580 


NGC 5548 


410 


Mrk 817 


330 


Mrk 509 


520 


NGC 7469 


360 



-51 - 



Table 4. Strength of Narrow-H/3 Component 



Object 


F(H/3 narrow )/F([Om] A5007) 


Ref. 1 


(1) 


(2) 


(3) 


PG 0026+129 


0.108 


1 


PG 0052+251 


0.120 


1 


Fairall 9 


0.121 


1 


PG 0953+414 


0.147 


1 


NGC 3227 


0.088 


2 


NGC 4151 


0.084 


2 


PG 1229+204 


0.080 


2 


PG 1307+085 


0.094 


1 


IC 4329A 


0.136 


1 


Mrk 279 


0.215 


1 


PG 1351+640 


0.330 


1 


PG 1411+442 


0.213 


1 


NGC 5548 


0.110 


2 


PG 1426+015 


0.120 


2 


PG 1613+658 


0.075 


1 


PG 1617+175 


0.192 


1 


PG 1704+608 


0.110 


2 


PG 2130+099 


0.306 


1 


NGC 7469 


0.181 


1 



References: 1: Marziani et al. 2002; 2: This work. 



Table 5. Time-Series Analysis 



Data Julian Dates F var F va , T r ccnt T poak 

Object Set (-2400000) r max (cont) (line) (days) (days) 

(1) (2) (3) (4) (5) (6) (7) (8) 



Mrk 335 


H/3 


49156-49338 







H/3 


49889-50118 





PG 0026+129 


Ha 


48836-51084 







H/3 


48545-51084 





PG 0052+251 


Ha 


48837-51084 







H/3 


48461-51084 







H 7 


48461-51084 





Fairall 9 


H/3 


50473-50665 







CivA1549 


50473-50713 







Lya 


50473-50713 





Mrk 590 


H/3 


48090-48323 







H/3 


48848-49048 







H/3 


49183-49338 







H/3 


49958-50122 





3C 120 


H/3 


47837-50388 





Akn 120 


H/3 


48148-48344 







H/3 


49980-50175 





Mrk 79 


H/3 


47838-48044 







H/3 


48193-48393 







H/3 


48905-49135 







H/3 


49996-50220 





PG 0804+761 


Ha 


48319-51085 







H/3 


48319-51085 







H 7 


48319-51085 





PG 0844+349 


Ha 


48319-51085 







H/3 


48319-51085 







H 7 


48319-51085 





Mrk 110 


H/3 


48953-49149 







H/3 


49751-49874 







H/3 


50010-50262 







HeiiA4686 


48953-49149 







He n A4686 


49751-49874 







HeiiA4686 


50010-50262 





PG 0953+414 


H/3 


48319-50997 







H 7 


48319-50997 





NGC 3227 


Ha 


47894-48045 







H/3 


47894-48049 







H/3 


48623-48776 






829 


± 





066 





066 





041 


17-3±1;1 


18l 5 6 


714 


± 





094 





049 





015 


12.8+2? 


131? 


796 


± 





047 





115 





092 


ii2.it22:? 


1081** 


814 


± 





047 





173 





079 


126.8l^j 


105±i 7 


646 


± 





063 





192 





101 


189.1±S;S 


1941H 


748 


± 





066 





199 





114 


103.7±i?i 




672 


± 





067 





199 





199 


94.21^-j 


75+ 7 2 


715 


± 





067 





328 





043 


18.2+tl 


19l| 


619 


± 





089 





385 





056 


: 3i.o±ig : « 


• 35+ 14 


898 


± 





026 





385 





115 


12.51°;° 


9+t 


709 


± 





099 





075 





063 


91 9+3.6 
^•^ — 2.8 


2ll| 


948 


± 





025 





103 





334 


14.4+J; 7 


17 +8 7 


874 


± 





257 





068 





104 


30.015° 


321* 


935 


± 





031 





153 





098 


29.51^ 


30 + 5 


646 


± 





070 





178 





095 


39.4l«i 


33+i° 


839 


± 





076 





039 





059 


48.6l? 2 6 8 


52l? 2 


940 


± 





025 





104 





112 


38.315° 


37+1 7 


684 


± 





108 





093 





057 


q 9+ 8 - 5 

y - — 8.0 


9lS 


870 


± 





053 





100 





081 


16.41?:? 


19111 


689 


± 





082 





095 





048 


16.41^ 


161J 


550 


± 





154 





077 





052 


• 5 7+ 2 - 7 

• °-'-4.3 


:8lf 


818 


± 





032 





176 





058 


195.3+^? 


1851^ 


831 


± 





031 





176 





061 


156.31™ 


120121 


685 


± 





070 





176 





080 


114.011?"? 


90+H 


705 


± 





068 





105 





067 


34.4±il| 


oc+8 
O0 -14 


759 


± 





055 





105 





096 


. o 9 +13.2 
• °-^-10.6 


= 6liS 


533 


± 





103 





105 





054 


or Q+ 74 - 6 


26l 7 7 


704 


± 





112 





102 





031 


25 1+ 5 ' 7 


26l?x 


829 


± 





506 





126 





062 


21.ll 6 ° 5 9 


21+ 7 2 


960 


± 





014 





323 





198 


34.4l£t 


25+5 


773 


± 





097 





102 





222 


: 0.712;? 


: 0.5l 2 2 i 


918 


± 





223 





126 





435 


: 2.4li 2 5 


: 1-5+21 


924 


± 





027 





323 





648 


: 8.5l^ 5 6 


: 0.012 ° 


666 


± 





101 





136 





054 


185.2l2?j 


2071 2! 


483 


± 





141 





136 





056 


188.91^1 


207! 25 , 


786 


± 





081 





128 





017 


19.0l?i 7 4 


201J1 


770 


± 





100 





151 





133 


5 5+ 14 ' 1 


4+14 
^-8 


599 


± 





106 





036 





075 


8.2±H 


n+6 

y -8 



Table 5 — Continued 



Data Julian Dates F var F var r ccnt r poak 

Object Set (-2400000) r max (cont) (line) (days) (days) 

(1) (2) (3) (4) (5) (6) (7) (8) 



NGC 3516 


Ha 


47894-48047 


0.755 ± 0.078 





289 





129 


13.31^ 


1311 




H/3 


47894-48047 


0.685 ± 0.086 





289 





110 


6.8i 6 3 -j 


7 +9 


NGC 3783 


H/3 


48607-48833 


0.666 ± 0.072 





192 





066 


10.4121 


9±i 




SiivA1400 


48611-48833 


0.457 + 0.081 





192 





108 


2.01?;? 


2.011;;; 




CivA1549 


48611-48833 


0.478 ± 0.074 





192 





069 


o o+l.O 
d -°-0.9 


4-011° 




HenA1640 


48611-48833 


0.582 + 0.068 





192 





154 


l-4l° j 


1 V. + 0.5 
1 -°-1.0 


NGC 4051 


H/3 


50183-50263 


0.763 ± 0.072 





059 





096 


5-8llJ 


6.0!?;! 




HenA4686 


50183-50263 


0.761 +0.086 





059 





114 


5.011;? 


5.oll;? 


NGC 4151 


H/3 


49305-49404 


0.837 ± 0.040 





057 





061 


3 1+1-8 
« 5 - 1 -1.3 


3.0ll;g 




Ha 


49305-49404 


0.815 + 0.038 





057 





033 


3.2+?'? 


i-5±|;g 




H/3 


47145-47360 


0.753 + 0.053 





058 





150 


n-6ll;l 


io.ol|° 




Ha 


47145-47360 


0.772 ± 0.054 





058 





096 


n-olll 


io.oll° 


PC 1211+143 


Ha 


48319-51000 


0.864 ± 0.038 





134 





103 


: KHU) • ::;;;, 


:65ll 3 




H/3 


48319-51000 


0.855 ± 0.047 





134 





121 


: 101.41&S 


: 29ll| 




H 7 


48319-51000 


0.832 ± 0.053 





134 





124 


: 139.411?;! 


: 133+lf 


PG 1226+023 


Ha 


48361-50997 


0.694 ± 0.071 





102 





077 


504.41 5 ™ 


458+ 2 ? 8 




H/3 


48361-50997 


0.787 ± 0.056 





102 





078 


otrtr n+79.4 
OOO.O_ 105 3 


337l^° 8 




H 7 


48361-50997 


0.726 ± 0.067 





102 





119 


304.4l«;g 


307ll 7 


PC 1229+204 


Ha 


48319-50997 


0.775 ± 0.050 





107 





081 


94.6lf 6 ;? 


5llf 




H/3 


48319-50997 


0.764 ± 0.055 





107 





125 


40.2lf 6 ; 3 


34+ 11 




H 7 


48986-50997 


0.731 ± 0.074 





096 





175 


11 6+ 211 

11.0_ 27 3 


7+28 
' -17 


NGC 4593 


Ha 


47894-48049 


0.776 ± 0.091 





114 





133 


3 2+ 5 - 6 

—4.1 


2+ 7 

z -2 




H/3 


47894-48049 


0.809 ± 0.099 





114 





185 


: 1.2l 9 5 ! 


• 2 +8 


PG 1307+085 


Ha 


49130-51000 


0.716 + 0.092 





118 





078 


164.8l? 32 5 9 


197lH 8 




H/3 


48319-51042 


0.725 ± 0.065 





113 





123 


121.9±&g 


iistS 




H 7 


48319-51042 


0.630 + 0.113 





113 





102 


192.0lf j 


222+^ 


IC 4329A 


H/3 


48643-48832 


0.616 + 0.115 





115 





084 


1 5+ 2 ' 7 


H? 


Mrk 279 


H/3 


50095-50289 


0.692 + 0.055 





092 





058 


17 n+4.0 
i '" 5 -4.0 


191^ 




H/3 


47205-47360 


0.731 +0.065 





071 





138 




• 12+ 12 

• lz -n 




Ha 


47205-47360 


0.758 + 0.055 





071 





067 


19 o+9.8 


1211? 


PG 1411+442 


Ha 


48319-51038 


0.830 ± 0.047 





105 





034 


103.2+^1 


7311? 




H/3 


48319-51038 


0.620 + 0.095 





105 





050 


135.4lf 7 J 


77 +161 
' ' -21 


NGC 5548 


H/3 


47509-47809 


0.813 + 0.030 





117 





091 


20.0li;| 


22 +2 




H/3 


47861-48179 


0.855 + 0.033 





129 





191 


18.9+24 


18+ 2 




H/3 


48225-48534 


0.659 ± 0.068 





090 





093 


16.H2:? 


16+| 




H/3 


48623-48898 


0.893 + 0.019 





168 





284 


11.21^ 


nil 




H/3 


48954-49255 


0.666 ± 0.044 





087 





057 


13.21^ 


1212 




H/3 


49309-49636 


0.792 ± 0.029 





104 





117 


i3.7±|:§ 


q+8 

9 -3 



Table 5 — Continued 



Data Julian Dates F var F var r ccnt r poak 

Object Set (-2400000) r max (cont) (line) (days) (days) 

(1) (2) (3) (4) (5) (6) (7) (8) 





H/3 


49679-50008 


0. 




H/3 


50044-50373 


0. 




H/3 


50434-50729 


0. 




H/3 


50775-51085 


0. 




H/3 


51142-51456 


0. 




H/3 


51517-51791 


0. 




H/3 


51878-52174 


0. 




HenA4686 


47509-47797 


0. 




SiivA1400 


47510-47745 


0. 




CivA1549 


47510-47745 


0. 




HeiiA1640 


47510-47745 


0. 




Cm] A1909 


47510-47745 


0. 




SiivA1400 


49097-49135 


0. 




CivA1549 


49060-49135 


0. 




HenA1640 


49097-49135 


0. 




Cm] A1909 


49097-49135 


0. 




H/3 


47212-47360 


0. 




Hq 


47212-47360 


0. 


PG 1426+015 


Ha 


48334-51042 


0. 




H/3 


48334-51042 


0. 


Mrk 817 


H/3 


49000-49212 


0. 




H/3 


49404-49528 


0. 




H/3 


49752-49924 


0. 


PG 1613+658 


H/3 


48397-51073 


0. 


PG 1617+175 


Hq 


48362-51085 


0. 




H/3 


48362-51085 


0. 


PG 1700+518 


H/3 


48378-51084 


0. 


3C 390.3 


H/3 


49718-50012 


0. 




HenA4686 


49718-50012 


0. 




Lya+Nv 


49718-50147 


0. 




CivA1549 


49718-50147 


0. 


Mrk 509 


H/3 


47653-50374 


0. 




He n A4686 


47653-50374 


0. 


PG 2130+099 


Ha 


48382-51084 


0. 




H/3 


48382-51084 


0. 




H 7 


48382-51084 


0. 


NGC 7469 


H/3 


50237-50295 


0. 




Hq 


50237-50295 


0. 



841 


± 





.032 


0. 


079 


0. 


071 


22.11-6 


23l 4 


.875 


± 





,019 


0. 


150 


0. 


.106 


16-711:1 


15+? 


753 


± 





,044 


0. 


105 


0. 


.125 


17.91?? 


17+ 4 


.800 


± 





,032 


0. 


.100 


0. 


097 


26.9+j ; 4 


2911 


802 


± 





,037 


0. 


149 





.058 


25.3«- 3 


251^ 


.845 


± 





,036 


0. 


166 


0. 


192 


6-611:1 


91* 


751 


± 





,053 


0. 


113 


0. 


208 


i4.5l 6 7 :° 


g+13 


.770 


± 





,051 


0. 


117 


0. 


052 


8.0l?;jj 


o+3 
8-3 


.550 


± 





,088 


0. 


321 


0. 


.185 


12.5+1' 5 


1311 


.737 


± 





,050 


0. 


321 


0. 


136 


io.ol?-° 


io±S 


752 


± 





,055 


0. 


321 


0. 


344 


3 9+ 1 ' 7 


2+t 


.571 


± 





,090 


0. 


321 


0. 


130 


27.9+|-| 


2611 1 


.814 


± 





,051 


0. 


141 


0. 


.131 


4-3±i:S 


4 +2 


893 


± 





,029 


0. 


141 


0.11 


6-811:2 


51? 


934 


± 





,020 





.141 


0. 


141 


i-9±8:i 


2.ol°j 


.880 


± 





,040 


0. 


141 


0. 


086 


14.111:1 


15tJ 


590 


± 





,090 


0. 


.034 


0. 


072 


8.011J 


7+5 
' -2 


.556 


± 





,083 


0. 


034 


0. 


031 


7-5±«:l 


6+^ 


.813 


± 





,084 


0. 


173 


0. 


070 


82.11^ 


58l 6 2 ° 


764 


± 





,102 


0. 


173 


0. 


084 


103.21^-j 


6811? 


.780 


± 





,075 


0. 


.135 


0. 


.051 


19-6+|-g 


21 + s 


.888 


± 





,059 


0. 


.096 


0. 


097 


15-sllJ 


17l| 


843 


± 





,068 





049 





054 


34.7+79 


36+?° 


.394 


± 





,095 





123 


0. 


060 


45.211?:? 


3411? 


756 


± 





,051 


0. 


191 


0. 


.072 


io4.8±|S:? 


94ll| 


701 


± 





,070 





.191 


0. 


108 


79.411?:° 


83lg 


.532 


± 





,095 


0. 


060 


0. 


021 


325.3lg;? 


326l^° 


.675 


± 





,063 


0. 


343 


0. 


072 


24.91?:? 


28l? 2 


424 


± 





,078 


0. 


343 


0. 


111 


oo '7+33.0 
Z °-'-25.3 


oo+34 
28_ 30 


777 


± 





,048 


0. 


343 


0. 


136 


6i.9iil:? 


62+ 32 
oz -36 


.871 


± 





,023 


0. 


.343 


0. 


294 


37.7+1?:} 


31111 


831 


± 





,025 





181 





.106 


82.3+11 


79+? 


466 


± 





,043 





.181 


0, 


336 


34.61? | 


15+?° 


.643 


± 





,040 


0. 


086 





.062 


210.9lH:l 


21211* 


.678 


± 





,042 


0. 


.086 





089 


168.6±?J : | 


156±H 


632 


± 





,071 


0. 


086 


0. 


082 


196.2+1?? 


220l« 


.654 


± 





,051 


0. 


150 


0. 


043 


4.7±g:S 


4-5lJ:° 


619 


± 





,074 


0. 


150 


0. 


018 


4-9111 


5.o±?-° 



Table 5 — Continued 



Data Julian Dates F var -FVar T cont r pca k 

Object Set (-2400000) r max (cont) (line) (days) (days) 

(1) (2) (3) (4) (5) (6) (7) (8) 



SiivA1400 50245-50293 
CivA1549 50245-50293 
HeiiA1640 50245-50293 



0.651 ±0.041 0.150 
0.669 ±0.035 0.150 
0.553 ±0.052 0.150 



0.189 1.8±g;l 2.0t° 5 5 
0.082 2.6±g;l 2.5±g;| 
0.135 0.6±°o.l 0.8±g;| 
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Table 6. Rest-Frame Lags, Line Widths, and Virial Products 



Tcent 7"pcak <^linc 

FWHM Cr ccn t ofine/ G 

Object Line (days) (days) (kms" 1 ) (kms" 1 ) (1O 6 M ) 

(1) (2) (3) (4) (5) (6) (7) 



Mrk 335 


H/3 


16.8±*"1 


18l 5 6 


917 ±52 


1629 ± 145 


2 76+ ' 85 

^- ' "J-0.76 




H/3 


12.5±|f 


13+? 


948 ± 113 


1375 ± 357 


2.2Bt\f 


PG0026+129 


Ha 


98.llij 


95lf 2 


1961 ± 135 


1117 ± 109 


73.7+l?' 7 




H/3 


m.ollti 


92lf 9 


1773 ± 285 


1719 ± 495 


68.211^ 


PG0052+251 


Ha 


163.7±||;1 


I68+77 


1913 ± 85 


: 2682 ± 453 


H7.0±«i 




H/3 


89.8l 24 4 ? 




1783 ± 86 


4165 ± 381 


55-711^ 




H 7 


8i.6±i»:S 


65l 10 


2230 ± 502 


5633 ± 3585 


79 9+ 40 ' 3 


Fairall 9 


H/3 


i7.4tij 


1812 


3787 ± 197 


6901 ± 707 


48.7±i°:S 




CivA1549 


: 29.6t£° 


. oq+13 
• 00 -19 


3201 ± 285 


4628 ± 1375 


: 59.3l 2 « : ° 




Lya 


ll-9±s.l 


9l 6 6 


4120 ± 308 


3503 ± 1474 


39.4l 2 °« 


Mrk 590 


H/3 


20.71^ 


20±1 


: 789 ± 74 


: 1675 ± 587 


. o C9+0.64 
• 0.58 




H/3 


14-0±1;| 


17+ 8 7 


1935 ± 52 


2566 ± 106 


10.25l« ; 21 




H/3 


29.21" 


3ll| 


1251 ± 72 


2115 ± 575 


8-9211: 81 




H/3 


28.8l^ 


291? 


1201 ± 130 


1979 ± 386 


s.iiil:?! 


3C 120 


H/3 


38.ll?J"l 


32+ 10 


1166 ±50 


2205 ± 185 


lO-lllI 


Akn 120 


H/3 


47.ll? 2 3 4 


50l? 2 


1959 ± 109 


5536 ± 297 


35.3_ 10 A 




H/3 


37.11" 


36+ 16 


1884 ± 48 


5284 ± 203 


25.71^ 


Mrk 79 


H/3 


9-0l 7 j 




2137 ±375 


5086 ± 1436 


8.05l 7 :? 2 




H/3 


16.lle1 




1683 ± 72 


4219 ± 262 


Q QQ+3-72 
o.8»_3 70 




H/3 


16-0±i;| 




1854 ± 72 


5251 ± 533 


10.76±|-|§ 




H/3 


: 5.6±i| 


:8l 12 


1883 ± 246 


2786 ± 390 


: 3.84l 2 :^ 


PG0804+761 


Ha 


183.61^ 


1741H 


2046 ± 138 


3155 ± 569 


i5o.o±H:| 




H/3 


146.9±HJ 


113115 


1971 ± 105 


2012 ± 845 


111 4+ 18 ' 6 




H 7 


107.ltg-° 


8511 


: 1256 ± 141 


: 2201 ± 295 


■■ 33.011?: 1 


PG0844+349 


Ha 


on q+13.7 
OZ.O_ 13 4 


oo+8 
00_ 13 


1625 ± 73 


2436 ± 329 


i6.7l 7: 2 




H/3 


: 3.0t»;3 




1448 ± 79 


2148 ± 612 


= 1-2111 




H 7 


27.HS:! 


24 + S 


2058 ± 218 


4946 ± 1085 


OO /l +58.2 

zz.^_ 31 9 


Mrk 110 


H/3 


24.3+83 


251^ 


1196 ± 141 


1494 ± 802 


6.78ll-il 




H/3 


20.4li d 3 5 


20l 7 2 


1115 ± 103 


1381 ± 528 


4.94+?- 7 ^ 




H/3 


oo o + li.9 
JJ -°-10.0 


24+3 


755 ± 29 


1521 ±59 


3.701J:? 8 




HeiiA4686 


:0.7tl:J 


: o.5tl-i 


2118 ± 138 


4975 ± 513 


• c iq+ 2 ' 72 




He n A4686 


• 2 5+ 4 ' 


: 1-4151 


2625 ± 85 


5253 ± 810 


• s 41 + 5 - 42 




HeiiA4686 


. s 9 +10.3 


: O.OlU 


2207 ± 83 


3416 ± 428 


■ 7 77+9-77 

■ '■' '-7.86 


PG0953+414 


H/3 


150.ll 2 2 ^ 


168±1| 


1306 ± 144 


3002 ± 398 


so.oli 3 : 1 




H 7 


153.1±l?;i 


168l 2 ° 


1299 ± 193 


2960 ± 589 


50.4115:2 


NGC 3227 


Ha 


18-9l!i r 3 


201J 2 


1977 ± 134 


3168 ±67 


14.43ll:|° 




H/3 


8-2+11 


n+6 

y -8 


1925 ± 124 


5138 ± 787 


5-9011:11 




H/3 




4+14 


2018 ± 174 


5278 ± 1117 


4.331^ 



Table 6 — Continued 



Tcent 7"pcak ^"linc 

FWHM CT ccnt af inc /G 



\j o j ecx 


Line 


(days) 


(days) 


(l m a — l~i 
^KHl S J 


[Km s j 


fin* \/f „ 1 


ft \ 




W 






W 


{') 


NGC 3516 


Ha 


13.2+IJ 


13+5 


2108 ±69 


4770 ± 893 


11.48+aif 




H/3 


6.7+5/* 

— d.es 


7 + s 


1837 ± 115 


3353 ± 310 


4.45+|-g 


NGC 3783 


CivA1549 


o + l.O 


4.0+1° 


2948 ± 160 


3691 ± 475 


6.3711:1? 




H/3 


10.2+11 


9+5 


1753 ± 141 


3093 ± 529 


6-141?;! 




HenA1640 


1.4+2-8 

— U . 


i.o+?-5 

— l.U 


3870 ± 162 


8008 ± 1268 


4.00+?'?r 

— 1 .OD 




SiivA1400 


2.0+?? 


2.0+!° 

— l.U 


3488 ± 161 


6343 ± 2021 


4-71+ 2 .! 2 . 


NGC 4051 


H/3 


5.8+?° 


5-5+11 


543 ± 52 


1072 ± 112 


0.34l°;l 6 2 




HenA4686 


5.0+1 ? 

— d. 1 


5.0+3° 

— O.O 


1690 ± 88 


3490 ± 918 


2.78+?li 

— 1.77 


NGC 4151 


Ha 


3.2+1?, 


1.5+1° 


2422 ± 79 


3156 ± 300 


71+2.23 
<->■ 1 -i 2 {){) 




H/3 


3.1+il 


3.0+|° 


1914 ±42 


4248 ± 516 


9 10+0.91 

Zi . 10 Q 90 




Ha 


11.0+|i 


10+I 


: 1721 ±47 


: 3724 ± 529 


6.361?" 




H/3 


H.5+I-; 


lots 


: 1958 ± 56 


: 5713 ± 6759 


8.63l|lt 


PG1211+143 


Ha 


: 92.5tS:S 


: 60+li 


: 2321 ± 231 


: 1425 ± 382 


q 7 n+37.6 

— 43. ( 




H/3 


: 93.8lg;f 


: 27 + \l 


1080 ± 102 


1317 ± 138 


21.4ll 1 4 




H7 


: 129.0+21 2 

-OO.O 


: 123+" 3 
—by 


1376 ± 157 


2014 ± 249 


47.7+lI'2 

— lo.U 


PG1226+023 


Ha 


435.5±«:S 


395+^° 


: 2075 ± 239 


: 1638 ± 424 


: 366l°t 




H/3 


306.8+2 8 ,- 5 . 
— yu.y 


291+2? 

— ol 


1777 ± 150 


2598 ± 299 


189+2? 

— 04 




H7 


262.8+g-| 


265+U 


1688 ± 142 


3274 ± 484 


146+H 


PG1229+204 


Ha 


89.0tg;i 


48+f 


1737 ± 118 


3229 ± 364 


52.4+21-3 

— /O.O 




H/3 


37.8l?£| 


32+1° 


1385 ±111 


3415 ± 320 


14.2+i° 2 6 




H7 


10.9+^-2 


— ID 


1540 ± 159 


3207 ± 875 


5-l+?i 2 i 


NGC 4593 


Ha 


3.2+5-f 


2 + 2 


1253 ± 90 


3399 ± 196 


0.97+1™ 




H/3 




:2+ 8 7 


1673 ± 169 


3769 ± 862 


0.66+5°i 

— z.oo 


PG1307+085 


Ha 


142.7ll? 5 2 


17lt?28 


1843 ± 98 


: 3084 ± 1041 


94.612I-5 




H/3 


105.6lf 6 ;° 


103l« 


1820 ± 122 


5058 ± 524 


68.312?;° 




H7 


166.3l« : ? 


192+1° 


1758 ± 193 


4278 ± 881 


100.3llf| 


IC 4329A 


H/3 


i-5±?:S 


H? 


2476 ± 226 


6431 ± 6247 


1.80+S?g 


Mrk 279 


H/3 


16.7±|;g 


18t$ 


1420 ± 96 


3385 ± 349 


6.6011™ 




Ha 


12 4+ 9 - 5 


12 +13 


1405 ± 266 


: 3408 ± 555 


4 70+4-08 
' a -4.94 




H/3 


. Q O + 20.9 

• 8 ' 2 -10.4 


: 1211? 


1900 ± 84 


: 5717 ±6815 : 


5.79ll 4 3 r 


PG1411+442 


Ha 


94.7t^;g 


6711? 


2437 ± 196 


1877 ± 375 


io9.8lg : * 




H/3 


124.3±|i;? 


Vi+148 
' J--19 


1607 ± 169 


2398 ± 353 


62.6l 3 3 ^ 


NGC 5548 


H/3 


19-7±1;| 


22+s 


1687 ± 56 


4044 ± 199 


io.9ll;l 




H/3 


I8.6+23 


18+ 2 


1882 ± 83 


4664 ± 324 


12 q +1 - 9 

—2.0 




H/3 


15.9 +2 2 ^ 


1612 


2075 ± 81 


5776 ± 237 


13-3+23 




H/3 




HIS 


2264 ± 88 


5691 ± 164 


11.012° 




H/3 


i3.o±i-2 


1211 


1909 ± 129 


: 2543 ± 605 


9-211:3 




H/3 




n+8 

y -3 


2895 ± 114 


7202 ± 392 


22.01^5 



Table 6 — Continued 



Tcent ''"peak ^"linc 

FWHM cT cant a'i inc /G 

Object Line (days) (days) (kms^ 1 ) (kms^ 1 ) (1O 6 M ) 

(1) (2) (3) (4) (5) (6) (7) 





H/3 


91 7+2.6 
zi - ' -2.6 


231* 


2247 ± 134 


6142 ± 289 


21-41^ 




H/3 


16.4+1? 


15+? 


2026 ± 68 


5706 ± 357 


i3.ili: 3 3 




H/3 


i7.5t; : g 


17+1 


1923 ± 62 


5541 ± 354 


12.61M 




H/3 


26.5t|:2 


2911 


1732 ± 76 


4596 ± 505 


15.5I?: 9 




H/3 


24.81^ 


251^ 


1980 ± 30 


6377 ± 147 


19.0l 2 j 




H/3 


6.51^ 


91S 


1969 ± 48 


5957 ± 224 


4.9+^ 




H/3 


14.31?:?, 


9l 3 3 


2173 ± 89 


6247 ± 343 


13.2l 5 6 j 




H/3 


7.8l2:| 


712 


3078 ± 197 


8047 ± 1268 


14 5+ 72 




HeiiA4686 




8l| 


2850 ± 151 


7338 ± 901 


12.4l|- 2 




SiivA1400 


12.3li: 4 


13l 2 


2576 ± 389 


6455 ± 3030 


15.912:2 




CivA1549 


9-8±i;l 


10l 3 


3842 ± 210 


6556 ± 878 


28.3±| : | 




HeiiA1640 


3-s±i:l 


2+ 4 


3897 ± 264 


9803 ± 1594 


ii-sl&S 




C m] A1909 


27.4l 5 5 j 


2611 1 


2360 ± 222 


5018 ± 1458 


29.8l| : ; 




SiivA1400 


4-3li:S 


41? 


4014 ± 253 


7044 ± 1849 


13.4l 3 j 




CivA1549 


6-71?:? 


51? 


3328 ± 104 


6868 ± 450 


14-51 2 : 2 




HenA1640 


1 Q+0-3 
1 - v -l).3 


2.ol°: 5 5 


4397 ± 154 


8929 ± 1571 


7 9+1.2 

' - z -1.3 




Cm] A 1909 


13.9111 


15l 2 


3227 ± 176 


4895 ± 1263 


28.2l 4 2 




Ha 




61? 


: 1694 ± 80 


: 3044 ± 381 


• 4 -I+ 3 - 2 


PG1426+015 


Ha 


75.5±%* 


53121 


4254 ± 290 


5450 ± 842 


267+ 114 
ZD ' -120 




H/3 


95.011?:? 


63l 9 2 ° 


3442 ± 308 


6323 ± 1295 


220l|5 


Mrk 817 


H/3 


19.0+ 3 ? 


20+5 


1392 ± 78 


3515 ± 393 


7 91+1.67 
' - zl -1.63 




H/3 


i5.3±l:J 


16l 3 4 


1971 ± 96 


4952 ± 537 


1 1 R9+ 3 ' 03 
ii -" JZ -2.90 




H/3 


33.6+?-^ 


35l?° 


1729 ± 158 


3752 ± 995 


19.63+5?3 


PG1613+658 


H/3 


40.lli 5 5 '° 2 


sol?! 


2547 ± 342 


7897 ± 1792 


50.7lll:| 


PG1617+175 


Ha 


94.2l«i 


84l2 4 6 


2483 ± 160 


3794 ± 780 


113.4+^? 




H/3 


7 1 r+29.6 
' L -°-33.7 




2626 ±211 


4718 ± 991 


96.2+f 7 J 


PG1700+518 


H/3 


251.81&? 


252lf 6 


1700 ± 123 


1846 ± 682 


142l 3 3 g 


3C 390.3 


CivA1549 


OJ - ' -14.6 


29l 2 ? 


5315 ± 165 1 


10357 ± 1494 1 


:196.7l« 4 : 2 




CivA1549 


or y+11.4 
od - ' -14.6 


29l 2 ? 


4400 ± 186 2 


8989 ± 2987 2 


134.81 44 :! 




H/3 


23.6l^'? 


27l?i 


3211 ±90* 


9630 ± 804 1 


■A7A+HI 




H/3 


23.61s: 2 . 


27l?x 


3105 ±81 2 


9958 ± 1046 2 


44.41 1 ! 2 :" 




He n A4686 


97 9+31.2 

z ' - z -24.8 


9 7 +32 
z ' -28 


2308 ± 155 1 


6046 ± 921 1 


■ OS Q+32.7 
• zo -°-26.0 




He ii A4686 


97 9+31.2 
Z ( .Z_ 24 g 


97+32 
z ' -28 


3305 ± 237 2 


8488 ± 1842 2 


58.olS: 5 




Lya 


58.6I272 


5911° 


4600 ± 141 1 


8225 ± 781 1 


•949 9+H5.2 
• z ^ z - z -113.2 




Lya 


58.6±1?:5 


5911° 


3952 ± 203 2 


8732 ± 985 2 


178.811U 


Mrk 509 


H/3 


79.61^ 


76+? 


1276 ± 28 


2715 ± 101 


25.3+ 2 '? 




HeiiA4686 


33.51?:? 


15l? 9 


2574 ± 85 


5864 ± 306 


43.31^° 


PG2130+099 


Ha 


198.4l2 2 j 


199!? 3 , 


1421 ± 80 


1574 ± 438 


78.2112:? 



-59- 



Table 6 — Continued 



Object 
(1) 


Line 
(2) 


Tccnt 

(days) 
(3) 


'peak 

(days) 
(4) 


Cline 

(kms -1 ) 
(5) 


FWHM 
(kms" 1 ) 
(6) 


CTcent ^"linc/ & 

(1O 6 M ) 
(7) 


NGC 7469 


H/3 

H 7 
CivA1549 

Ha 

H/3 
HenA1640 
SiivA1400 


i58.it":? 

184.5^'j 

9 c+0.3 
^■°-0.2 

4-711:1 

4.5l° i 

0- 6l° j 

1 y+0.3 

1 - '-0.3 


1471™ 
207±|| 

9 4 +0.5 
Z -*-0.5 

4 7+1-9 

4.3l° -j 
0.8l° j 
l-9l°o-j 


1623 ± 86 
1836 ± 191 
2619 ±118 

1164 ±68 
1456 ± 207 
3723 ±113 
3495 ± 269 


2912 ± 231 
2661 ± 481 
4305 ± 422 
1615 ±119 
2169 ± 459 
10725 ± 1697 
6033 ± 1112 


81 Q+17-6 

121 4+ 401 

q qc+0.47 
O.OO_ 45 

1 94+0.44 
0.36 

1 67+0- 89 
4 17+0.94 

-0.92 



1 Measurements for blended line. 

2 Measurements based on unblended half of line. 
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Table 7. Virial Relationship Fits 



Free Slope Fixed Slope b = —0.5 

Data Set a b \l a xl 



(1) (2) (3) (4) (5) (6) 



NGC 5548: RMS (Fig. 3) 



°"linc vs - feent 


3. 


926 ± 0. 


082 


-0. 


.497 ± 0. 


.070 


4.04 


3 


.930 ± 


.016 


3.86 


FWHM VS. Tccnt 


4. 


262 ± 0. 


133 


-0. 


427 ± 0. 


112 


4.25 


4 


.346 ± 


021 


4.12 


CT line ''"peak 


4. 


.033 ± 0. 


138 


-0. 


.583 ± 0. 


117 


2.96 


3 


.936 ±0. 


.018 


2.88 


FWHM vs. r peak 


4 


187 ±0. 


102 


-0. 


378 ± 0. 


086 


2.09 


4 


,328 ± 


.020 


2.13 



NGC 5548: Mean 1 (Fig. 4) 



°"linc vs - feent 


4.058 ±0.135 


-0.545 ± 0.114 


7.60 


4.005 ±0.017 


7.30 


FWHM VS. Tccnt 


6.202 ± 1.602 


-2.277 ± 1.357 


14.42 


4.108 ±0.041 


29.15 


CT line ''"peak 


4.264 ±0.268 


-0.721 ± 0.228 


4.25 


4.005 ± 0.020 


4.44 


FWHM VS. T p c ak 


8.035 ±3.761 


-3.882 ± 3.209 


5.95 


4.081 ±0.046 


19.23 


NGC 5548: Mean 2 (Fig. 5) 












""line vs - Tcent 


3.917 ± 0.072 


-0.422 ± 0.060 


4.70 


4.009 ± 0.014 


4.74 


FWHM VS. rccnt 3 


3.870 ± 0.066 


-0.088 ± 0.060 


70.25 


4.341 ± 0.033 


54.68 


°"linc vs - ''"peak 


3.987 ±0.112 


-0.483 ± 0.095 


2.78 


4.006 ±0.015 


2.65 


FWHM vs. r peak 3 


3.871 ±0.057 


-0.092 ± 0.053 


46.90 


4.328 ± 0.034 


31.11 


NGC 3783: RMS (Fig. 6) 












"line vs - feent 


3.690 ± 0.042 


-0.428 ± 0.063 


0.28 


3.735 ±0.021 


0.28 


FWHM VS. -r ccnt 


3.925 ±0.098 


-0.497 ±0.144 


0.74 


3.927 ±0.037 


0.50 


0"linc vs. r pcak 


3.696 ± 0.078 


-0.446 ± 0.122 


0.52 


3.727 ±0.030 


0.37 


FWHM vs. r peak 


3.894 ± 0.069 


-0.473 ± 0.104 


0.25 


3.910 ± 0.026 


0.17 


NGC 7469: RMS (Fig. 7) 












""line vs - feent 


3.786 ±0.114 


-0.952 ± 0.248 


1.19 


3.592 ±0.042 


3.36 


FWHM vs. r cont 


4.066 ± 0.093 


-1.119 ± 0.195 


0.50 


3.772 ±0.067 


3.94 


"line vs - fpeak 


3.737 ±0.105 


-0.918 ± 0.209 


0.69 


3.553 ±0.049 


2.32 


FWHM vs. r peak 


4.028 ± 0.084 


-1.113 ± 0.164 


0.28 


3.742 ± 0.073 


3.36 


3C 390.3: RMS 4 (Fig. 8) 












CT linc vs - r ccnt 


2.083 ± 1.027 


±1.009 ± 0.686 


0.72 


4.339 ±0.097 


5.09 


FWHM VS. -r cont 


4.818 ±0.991 


-0.566 ± 0.659 


1.20 


4.718 ±0.045 


0.80 


3C 390.3: RMS 5 (Fig. 8) 












"line vs - feent 


2.660 ± 0.547 


±0.598 ± 0.365 


0.64 


4.300 ± 0.076 


3.10 


FWHM VS. -r cont 


4.197 ±0.108 


-0.149 ± 0.070 


0.13 


4.723 ± 0.036 


0.39 



includes narrow-line components. 
2 Narrow component of H/3 removed. 
3 Free-slope fit is unweighted. 
4 Full line width. 
5 Unblended wing of line. 
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Table 8. Adopted Virial Products and Derived Black-Hole Masses 



(cTcentffine/G) Mbh 1 log \L X (optical) 
Object (1O 6 M ) (1O 6 M ) (ergs s" 1 ) 



I 1 ) 


(+J 






(4) 






Mrk 335 


2.58 ± 0.67 


14.2 ± 3.7 


43. 


.00 


1 

zb 


o 


.04 


PG 0026+129 


71.4 ± 17.4 


393 ± 96 


45. 


no 
.Uz 


_|_ 


o 


,06 


PG 0052+251 


67.1 ± 13.8 


369 ± 76 


44. 


yo 


1 

± 





.08 


Fairall 9 


46.3 ± 10.1 


255 ± 56 


41. 


.zo 


_!_ 

□Z 





.05 


Mrk 590 


8.64+1.34 


47.5 ± 7.4 


43. 




l 

± 





.09 


3C 120 


10 1+iT 

—4. 1 


55 5+lli 

< -" J -^'— 22.5 


44 


1 *7 


1 

± 


o 


.08 


Akn 120 


27.2 ± 3.5 


150 ± 19 


44. 


O Q 
ZO 


_J_ 

zt 





.09 


Mrk 79 


9.52 ± 2.61 


52.4+ 14.4 


43. 


70 
. / Z 


_|_ 





,05 


PG 0804+761 


126 ± 15 


693 ± 83 


44. 


n 1 
y4 


1 

± 





,08 


PG 0844+349 


16.80 ± 7.06 


92.4 ± 38.1 


44. 


.00 


1 

zb 





.04 


Mrk 110 


4.57 ± 1.1 


25.1 ± 6.1 


43. 


*70 
. / Z 


i 

zb 





.09 


PG 0953+414 


50.1 ± 10.7 


276 ± 59 


45. 


.zz 


_i_ 
zb 





,06 


NGC 3227 


7.67 + 3.90 


42.2 ± 21.4 


42. 


QQ 

oo 


i 

zb 





.04 


NGC 3516 


7.76 ± 2.65 


42.7 ± 14.6 


42. 


8.8. 


_|_ 
zt 





.13 


NGC 3783 


5.42 ± 0.99 


29.8 + 5.4 


43. 


Of! 
.ZD 


_i_ 
zb 





.04 


NGC 4051 


0.348 ± 0.142 


1.91 ± 0.78 


41 


yo 


i 

zb 





.03 


NGC 4151 


2.42 ± 0.83 


13.3 ± 4.6 


42. 


.oo 


zb 





.23 


PG 1211+143 


: 26.6 + 8.0 


:146±44 


44. 


.75 


± 





,07 


PG 1226+023 


161+34 


886 ± 187 


45 


.96 


± 





,05 


PG 1229+204 


13.3 + 6.4 


73.2 ± 35.2 


44 


.08 


± 





.05 


NGC 4593 


u.»/o_ 1 264 


5 36+ 9 ' 37 


43. 


.09 


± 





.14 


PG 1307+085 


80.0 ± 22.3 


440 ± 123 


44. 


.88 


± 





.04 


IC 4329A 


: 1.80±1;;| 


■q qn +17 - 88 


43 


.32 


± 





.05 


Mrk 279 


6.35 + 1.67 


34.9 + 9.2 


43. 


.88 


± 





.05 


PG 1411+442 


80.5 + 26.5 


443 ± 146 


44. 


.63 


± 





,04 


NGC 5548 


12.20 + 0.47 


67.1 + 2.6 


43. 


.51 


± 





,11 


PG 1426+015 


236 ± 70 


1298 ± 385 


44. 


.72 


± 





,07 


Mrk 817 


8.98 ± 1.40 


49.4 + 7.7 


43. 


.82 


± 





.05 


PG 1613+658 




279 ± 129 


44. 


.98 


± 





.05 


PG 1617+175 


108 ± 25 


594 ± 138 


44 


.48 


± 





.08 


PG 1700+518 


1 49+ 33 


781+lH 


45 


.63 


± 





.03 


3C 390.3 


52.2 ± 11.7 


287 ± 64 


43. 


.95 


± 





.07 


Mrk 509 


26.0 ± 2.1 


143 ± 12 


44. 


.28 


± 





.08 


PG 2130+099 


83.0 + 10.0 


457 ± 55 


44. 


.46 


± 





.04 


NGC 7469 


2.21 ± 0.25 


12.2 + 1.4 


43. 


.72 


± 





.02 



1 Assuming / = 5.5 



